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Abstract 


Circumstellar disks of planetary debris are now known or suspected to closely orbit hundreds 
of white dwarf stars. To date, both data and theory support disks that are entirely contained within 
the preceding giant stellar radii, and hence must have been produced during the white dwarf phase. 
This picture is strengthened by the signature of material falling onto the pristine stellar surfaces; 
disks are always detected together with atmospheric heavy elements. The physical link between 
this debris and the white dwarf host abundances enables unique insight into the bulk chemistry 
of extrasolar planetary systems via their remnants. This review summarizes the body of evidence 
supporting dynamically active planetary systems at a large fraction of all white dwarfs, the rem- 
nants of first generation, main-sequence planetary systems, and hence provide insight into initial 
conditions as well as long-term dynamics and evolution. 
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1. Introduction 


A circumstellar disk or ring is a relatively planar structure containing particulate matter and or- 
biting a star. Disks are distinct from spherical clouds or envelopes of dust (and gas) that typically 
surround protostellar objects and some giant stars. Circumstellar disks appear at every stage of 
stellar evolution, though the origin of the orbiting material is not always clearly understood. Pre- 
main sequence stars accrete material from a disk that is the flattened remnant of the cloud from 
which they formed. Disks found at young stars in subsequent evolutionary stages are the likely 
sites of planet formation, migration, and sometimes destruction. Mature main-sequence stars ex- 
hibit dusty disks owing to recent or ongoing energetic collisions among orbiting asteroid or comet 
analogs — sometimes referred to as the “Vega phenomenon” after the first star observed to have or- 
biting, non-stellar material. Several first-ascent and asymptotic giant stars are also known to have 
circumstellar disks, though their origins are still debated with hypotheses ranging from debris in a 
cold cometary cloud to consumed stellar companions. It can be useful to consider debris disks as 
the entire planetary system, minus the planets; for detailed observational and theoretical reviews 


of debris disk phenomena, see|Zuckerman] (2001) and {Wyatt} (2008). 
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The Solar System has a debris disk, part of which can be observed as zodiacal light, and has 
two major components in the asteroid and Kuiper belts, which provide a basic model for extrasolar 


debris disks. Beginning with the discovery of circumstellar dust at Vega 1984), 


it has gradually emerged that the majority of dusty main-sequence stars have (enhanced) disks 
analogous to the relatively icy and cold Kuiper belt 2006). In 
contrast to the zodiacal cloud of the Solar System, detectable extrasolar debris disks tend to be 
relatively clear of dust in their inner regions, with a few notable exceptions where warm analogs to 


the main asteroid belt have been detected (Backman et al.|{2009} 2013), and hot, transient 
dust likely due to recent giant impacts in the terrestrial zone (Song et al.|/2005; 2010b; 
2014), 


Disks orbiting white dwarfs are a recently recognized phenomenon, with nearly 40 identifi- 
cations since 2005, where the bulk required discovery and characterization from space. These 
evolved systems are a sensitive probe of the frequency and long-term dynamics of planetary sys- 
tems formed around intermediate-mass stars, and provide a unique window on planetesimal chem- 
istry. A journal review of this field is compelling as none exist outside of book chapters, and there 
is now more than a decade of infrared and optical observations of disks, including sensitive space- 
based data provided by the Spitzer Space Telescope. Published theoretical studies of disk formation 
and evolution have increased in the last several years, and an overall comparison of theory with 
observation is needed. 

It is noteworthy that the known circumstellar disks orbiting white dwarfs have no analog among 
main-sequence or giant stars. However, as described in this review, it is likely that white dwarf 
planetary systems resemble, at least in part, both solar and extrasolar planetary systems, including 
their planetesimal belts, while the disks themselves most resemble the planetary rings of the Solar 
System. This article focuses on the debris disk phenomenon for white dwarfs from a primarily 
observational perspective, but also discusses the relevant theoretical models; for an in-depth review 
of post-main sequence, planetary system evolution theory, see (2016). In the context of 
planetary system assembly and origins, it is fundamental that scientists can (indirectly) obtain 
extrasolar planetesimal compositions from the heavy element abundances of debris raining onto 
white dwarf atmospheres; for a review of this topic, see Jura & Young] (2014). Partly in parallel, 
the aim of this article is to review the physical data and theory (from the stellar surface outwards) 
that makes the above scientific analysis possible, and to provide additional and broader context for 
white dwarf planetary systems. 

The review begins in §2 with a historical perspective, placing white dwarfs into the context 
of debris disk observations and introducing the pollution phenomenon. Next, §3 expounds on 
the singular influence of Spitzer and the rapid gain of observational knowledge during this era. 
§4 discusses the standard accretion processes and diagnostics, summarizes the overall observed 
properties of known polluted white dwarfs, and reviews the various surveys that have provided the 
best statistical inferences. Some possible physical origins and evolution of the disks themselves 
are discussed in §5, and a summary is provided in §6. 


2. History and Background 


Owing to the bright and variable infrared sky, astronomical observations at these wavelengths 
where dust emission is most prominent are challenging; the historical detection of a debris disk 
orbiting the 5" brightest star required the Infrared Astronomical Satellite (IRAS). Furthermore, 
because any (heated) circumstellar material derives its luminosity from the central star, the small 
radii of white dwarfs implies orbiting dust can be remarkably hard to identify, as they are typically 
~ 10+ times less luminous than their main-sequence counterparts. Despite the unprecedented 
sensitivity and all-sky coverage of JRAS, not a single white dwarf was detected at any wavelength 


(Shipman\|1986). This is perhaps unsurprising given that JRAS was only sensitive to point sources 
brighter than 500 mJy at its shortest wavelength bandpass of 12ym, while the brightest white dwarf 


in the sky (Sirius B) emits around 5 mJy there (Skemer & Close} |2011). 

Nevertheless, the first dusty white dwarf was discovered within four years of the Vega dis- 
covery. In strange contrast, however, white dwarfs are a relatively recent addition to the classes 
of stellar objects known to host circumstellar disks because additional examples remained hidden 
until 2005. It is noteworthy that if white dwarfs only hosted debris analogous to main-sequence 
stars, their disks would still remain largely undetected today, even in the presence of dedicated and 
relatively advanced, infrared space missions. Circumstellar dust at white dwarfs has not only been 
historically hard to detect, it remains considerably difficult or impossible to study in detail. Only 
future space-based, and innovative ground-based, facilities can significantly impact the status quo 
for white dwarfs in the infrared. 

Of particular relevance for white dwarfs is that infrared excess emission associated with a star 
can arise from irradiated circumstellar dust or from a self-luminous companion. Owing to their 
compact nature, white dwarfs can be easily outshone in the near-infrared by low mass stellar and 
even brown dwarf companions. Recognizing this fact, R. Probst searched a large sample of nearby 
white dwarfs for photometric excess at 1 — 2,um to study the luminosity function of the lowest 
mass stars and brown dwarfs, as companions (1983). This was the first search for 
infrared excess at white dwarfs, and Probst deserves the credit for an insight that is now taken for 
granted, and which fostered an abundance of subsequent infrared work on white dwarfs, including 
the first discovery of circumstellar dust}?| 


2.1. The Discovery of Infrared Excess at G29-38 
Inspired by the work of Probst and armed with a rapidly evolving set of detectors at the NASA 


Infrared Telescope Facility (IRTF), identified the first white dwarf 
with infrared excess emission that was not associated with a stellar companion; G29-38. Photo- 
metric observations at three bandpasses longward of 2 um revealed T ~ 1200 K emission in excess 
of that expected for the stellar photosphere of this relatively cool white dwarf (Figure (I). 
However, as the observations were driven by a search for substellar companions, the infrared 
emission was initially attributed to a spatially unresolved brown dwarf. Heated circumstellar dust 


*For historical completeness and accuracy, |Probst|(1983) detected K-band emission from the disk orbiting G29-38 
sometime between 1978 and 1980. Although the single bandpass data were insufficient to confidently identify the 
excess emission, he gathered the first photons from extrasolar planetary material(!), several years before the Vega and 
B Pictoris disks were detected by JRAS. 
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Figure 1: The prototype dusty white dwarf G29-38. Upper panel: The discovery of infrared excess made at the IRTF. 
The plot is a reproduction of the figure presented in [Zuckerman & Becklin| (1987) and shows their measurements 
with errors in five infrared bandpasses labelled in the figure J (1.25 ym), H (1.65 ym), K (2.2 um), L (3.5 um), and M 
(4.8 um). The dotted line is a blackbody consistent with the Te¢ = 11500 K stellar photosphere in the infrared. Lower 
panel: The detection of photospheric calcium in the optical (Koester et al.|[1997), at the same atomic transition that 
is strongest in the Sun (the Fraunhofer K line). The thin line is the actual data while the thick line is a model fit. The 
Ca Il absorption appears rather weakly beside the pressure-broadened Balmer lines in this hydrogen-rich star. 


was not favored due to the likelihood of rapid dissipation by radiation drag, with 
prophetically remarking that if a debris disk were orbiting sufficiently close to 
attain such high temperatures, then spectral signatures of accretion should be seen. 

Over the next few years, the infrared emission at G29-38 was studied intensely by many 
groups, and its unique properties sparked interest across several subfields of astrophysics research: 
brown dwarf and planet hunters, asteroseismologists, infrared astronomers, and white dwarf pun- 
dits. Observational evidence gradually began to disfavor a brown dwarf as the source of infrared 
emission. First, some of the first near-infrared imaging arrays revealed G29-38 to be a point 
source in several bandpasses. Second, near-infrared spectroscopy measured a continuum flux 


source (‘Tokunaga et al.,|1988), whereas a very cool atmosphere was expected to exhibit absorp- 


tion features. Third, the detection of optical stellar pulsations echoed in the near-infrared were 
difficult to reconcile with a brown dwarf secondary (1990). 
Fourth, significant 10 um emission was detected at G29-38, a few times greater than expected for a 
cool object with the radius of Jupiter, essentially ruling out the brown dwarf companion hypothesis 


(Tokunaga et al.|/1990}/Telesco et al.}/1990). 


A decade after the discovery of its infrared excess, the optical and ultraviolet spectroscopic 


detection of multiple metal species in the atmosphere of G29-38 1997) made it 


clear the star is currently accreting from its circumstellar environs (Figure[I). 


2.2. White Dwarf Atmospheric Pollution 


It would not be possible to accurately review circumstellar disks at white dwarfs, including the 
pivotal role of G29-38, without introducing the phenomenon of atmospheric metal contamination. 
Rather than being a central focus of this review, the topic of photospheric metals will be discussed 
in the context of, and its relation to circumstellar material around white dwarfs. 

The origin and abundances of photospheric metals in isolated white dwarfs has been an astro- 
physical curiosity dating back to the era when the first few white dwarfs were finally understood 


to be subluminous via the combination of spectra and parallax 1919). In a half page 
journal entry, (1917) reported (via W. Adams) that his accidentally discovered faint 
star with large proper motion had a spectral type of “about FO”, almost certainly based on its strong 


Call H and K absorption features (Figure|2) 
Any primordial heavy elements in white dwarfs can only be sustained in their photospheres 
for the brief period while the star is still rather hot and contracting, and then only to a certain 


degree (Chayer et al. |1995} 2014). For Tepp < 25,000 K, gravitational settling 


is enhanced by the onset of convection and atmospheric heavy elements sink rapidly in the high 
surface gravities of white dwarfs leaving behind 
only hydrogen or helium. Downward diffusion timescales for heavy elements in cool white dwarfs 
are always orders of magnitude shorter than their evolutionary (cooling) timescales 
jal.||1986), and thus external sources must be responsible for any photospheric metals. 


3With hindsight, the detection of metals in vMa 2 in 1917 was the first observational evidence of planetary systems 


beyond the Sun (Zuckerman| /2015). 


2.3. Stellar, Interstellar, or Circumstellar Matter 


There are only three possible sources for the atmospheric metals seen in isolated white dwarfs: 
1) stellar material either primordial or fallback, 2) accretion from the interstellar medium, or 3) 
infall from an intrinsic circumstellar environment (i.e. a remnant planetary system). Regardless 
of origin, the deposition of heavy elements necessary to enrich a white dwarf atmosphere may be 
facilitated by the formation of a disk, and thus photospheric metals and circumstellar disks are 
likely to have a fundamental physical connection. 

The argument against a primordial origin for atmospheric heavy elements is the fact that they 
rapidly sink in the high gravity environment of white dwarfs. This theoretical expectation is cor- 
roborated by high-resolution spectra that typically exhibit pure hydrogen or helium compositions 


(Zuckerman et al.| {2003} 2005). Once the stellar ember has cooled below 25 000 K, 


metals can no longer be radiatively supported and diffuse from the atmosphere in a matter of 
days to years (2009). At least 20 Myr is required to reach this temperature, and rough] 
0.5 Myr is spent as a hot and luminous star with Tog ~ 10°K, L~50 Los 
During this time, any circumstellar gas and small particles (a < 1 cm) within several AU, such as 
dust formed within a stellar outflow during the giant branches, will be effectively removed from 
the system by engulfment, sublimation, radiation pressure, or drag forces. Because white dwarf 
disks and pollution are observed at 30 — 600 Myr cooling ages, the circumstellar material must 
derive from relatively parent large bodies, or originate beyond several AU. 

Historically, accretion from the interstellar medium was the most widely accepted hypothesis 


for metals observed in white dwarfs. This interpretation was at least partly due to the fact that until 
1997, all but one case of photospheric metals had been detected in relatively cool, helium-rich stars 


(spectral type DZ or DBZ; [Sion et al.]1990b). Such white dwarfs have relatively transparent outer 
layers that aid the detection of heavy elements (compare e.g. Figures|1}and|2), and relatively deep 
convection zones with metal diffusion timescales up to 10° yr The relatively long, 
atmospheric residence times in these stars allow for the possibility that their metals could be the 
residuals of interstellar cloud encounters [1993alb). 

The interstellar accretion model faced two major obstacles even before dust disks emerged as 
a class of objects associated with metal-rich white dwarfs. Fundamentally, the nearly universal 
lack of hydrogen in DBZ stars has been a glaring problem, even for the larger class of helium 
atmosphere white dwarfs (1993). In essence, 
if white dwarfs accrete sufficient material of roughly solar composition, they should all (even- 
tually) exhibit hydrogen-dominated outer layers, whereas a significant fraction of all cool white 
dwarfs are hydrogen-deficient stars, including many with metals. Another serious challenge came 
from the detection of metals in stars with hydrogen-rich atmospheres (spectral type DAZ), where 
the heavy element diffusion timescales are only days to years. The first was confirmed in 1983 


(Lacombe et al.}/1983), but remained anomalous until 1997 when two additional examples were 
found, including G29-38 (Holberg et al.| {1997} 1997). A targeted search brought 
further detections the following year (Zuckerman & Reid\ |1998), and because these stars must 


have experienced recent or ongoing accretion, the DAZ white dwarfs implied stringent constraints 
for interstellar accretion scenarios. 
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Figure 2: The prototype metal-enriched white dwarf vMa 2. Upper panels: The 1917 photographic plate spectrum of 
vMa 2, and plate sleeve with handwritten notes by observer W. S. Adams (van Maanen]|1917); these historical records 
were kindly supplied by Carnegie Observatories, who maintains the Mt. Wilson archives (J. S. Mulchaey 2015, private 
communication). Box highlighted are the strong Ca II H and K Fraunhofer lines, and are relatively easy to see in the 
century old spectrum. In contrast to hydrogen-rich atmospheres (e.g. Figure (Ip, metal absorption features in helium- 
rich stars can be prominent, and often dominate their optical spectra. Lower panel: An unpublished optical spectrum 


of vMa2 taken with UVES on the Very Large Telescope (VLT). All salient features are absorption due to Fe, Mg, or 
Ca. 


These observational developments and outstanding theoretical issues led some to consider cir- 
cumstellar accretion models. The earliest of these was the cometary impact model put forth by C. 
Alcock, and represented the first recognition of white dwarf pollution as a potential exoplanetary 


system signature (Alcock et al.\|1986). A clear strength of cometary impacts is greater consistency 
with the hydrogen-poor nature of the accreted material in numerous DBZ stars, although it was 


soon noted that volatile-depleted planetary material was more appropriate [1990a). 
Other difficulties for the comet cloud model came from two DAZ stars; G238-44 has a sufficiently 
short metal diffusion time that requires a continuous rain of comets (1997), while 
the circumstellar debris around G29-38 could not be generated by an impact 
1998), 


2.4. The Asteroid Disruption Model 


Two major developments occurred in 2003, representing a historic crossroads for white dwarf 


disks and pollution just prior to the launch of the Spitzer Space Telescope 2004). 
(2003) published a high-resolution spectroscopic study of photospheric metals 


in a large sample of cool, hydrogen-rich white dwarfs, finding several new examples of DAZ 
stars and establishing a minimum frequency of 20% for the pollution phenomenon in single stars. 
The results were shown to be inconsistent with the cometary impact model in the following way: 
it was predicted that there should be sharp rise in the number of DAZ with decreasing calcium 
abundance, owing to more frequent impacts of smaller objects (1986), but instead 
a relatively flat distribution was found in calcium abundance over several orders of magnitude. 
Neither could the DAZ sample be reconciled with accretion from interstellar medium, due to the 
lack of dense clouds in the Local Bubble, a result also found previously for a modest number of 
DBZ stars (Aannestad et al.|(1993). 

With the sole exception of G29-38 there was a distinct lack of reliable (infrared) data on the 
circumstellar environments of white dwarfs, yet a growing profusion of stars contaminated by 
metals, and problems with existing hypotheses. In a short but seminal paper, Jura|(2003) modeled 
the observed properties of G29-38 by invoking a tidally-destroyed minor planet (i.e. asteroid) that 
generates an opaque, flat disk of dust analogous to the rings of Saturn. Rather than impacting 
the star, an asteroid perturbed into a highly eccentric orbit passes within the stellar Roche limit 
and is torn apart by gravitational tides (Debes & Sigurdsson} 2002). Ensuing collisions reduce the 
fragments to rubble and dust, dynamically relaxing into a flat configuration over multiple orbital 
timescales. The closely orbiting dust is heated by the star producing an infrared excess, and 
slowly rains down onto the stellar surface, polluting its otherwise-pristine atmosphere with heavy 
elements. 

The tidally disrupted asteroid model, and accretion from a remnant exoplanetary system, is 
now the standard framework for white dwarf metal pollution. Interestingly and years prior, 
deduced that the closely orbiting dust at G29-38 must have resulted from a 
relatively recent, catastrophic event near the white dwarf, such as a near collision with an asteroid 
or comet. This qualitative model for the infrared emission at G29-38 invoked cm-sized pebbles as 
these would not be rapidly destroyed near the star, but the flat disk model achieves this by shielding 
the bulk of material from the full stellar radiation field. In the sections that follow, it will be shown 


that this model has been strongly corroborated over the past decade by myriad infrared, optical, 
and ultraviolet data. 


3. The Spitzer Revolution and Era 


Fortunately, the last decade has seen a profusion of white dwarf disk discoveries, due primarily 
to the 2003 launch and unprecedented infrared performance of Spitzer. Important roles have also 
been played by projects with large sky coverage such as the Sloan Digital Sky Survey (SDSS; 
Gunn et al.[2006), the Two Micron All Sky Survey (2MASS;|Skrutskie et al.[2006), and the Wide 


Field Infrared Explorer (WISE; \Wright et al.|/2010). The number of confirmed metal-rich white 
dwarfs with circumstellar dust consisted of a single system from 1987 until 2005, but has grown 


to nearly three dozen (and counting) within a decade. 

Although there have been some infrared disk detections from the ground longward of 2.5 um 
(Tokunaga et al, [1990 Becklin et al, 2005} Melis et al 20TT, most systems are (far) too faint to 
be studied in this manner, and all confirmed dusty white dwarfs at the time of writing have been 
observed with Spitzer. As demonstrated in this section, it is only beyond 2.5 um that the bulk of 
infrared excesses become salient, and thus the final NASA Great Observatory played a leading and 
fundamental role in this field, and specifically in connecting atmospheric metals in white dwarfs 
to circumstellar debris. 


3.1. Basic Dust Properties 


Within the first two cycles of Spitzer observations, the only two known white dwarfs with cir- 
cumstellar dust (G29-38 and GD 3624) were observed with all three cryogenic instruments 


(2007b). The Infrared Array Camera (IRAC; obtained 
photometry at 3.6, 4.5, 5.7, and 7.9 um, the Infrared Spectrograph (IRS; ac- 
quired R ~ 100 spectra between 5 and 15 um, and the Multi-band Imaging Photometer for Spitzer 
(MIPS; Rieke et al.[2004) measured fluxes at 24 um for the two earliest examples shown in Figure 
and later for a number of similar white dwarf targets prior to the loss of cryogen in 2009. 
There are a few remarkable, shared properties that can be seen in these two infrared datasets, 
two of which turn out to be ubiquitous (to date) among dusty and polluted white dwarfs. First are 


the strong silicate emission features (Reach et al.|/2005a} 2007b). In fact, the emission 


observed at GD 362 is among the strongest silicate features seen towards any mature star (compare 


b.2. 2015} 2010b). Second, there is a notable lack of relatively cool dust, 


unlike the bulk of all debris disks orbiting main-sequence stars that exhibit infrared excess only at 
24 or 70 um MIPS observations (2008). While models indicate dusty stars 
such as Vega harbor substantial planetary material consistent with orbital regions analogous to the 
Kuiper Belt, those same (blackbody) models would place 1000 K dust orbiting these two white 
dwarfs within 0.7 Re (roughly 50 stellar radii). The next wave of disks discovered and character- 
ized photometrically by Spitzer strongly suggested that a lack of cool material was a hallmark of 


‘Tronically, the first dusty white dwarf to be discovered in the Spitzer era, and the second known example of the 


phenomenon, was identified from the ground (Becklin et al.||2005 2005). 
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Figure 3: The first two dusty and metal-enriched white dwarfs observed with all three cryogenic instruments aboard 


Spitzer (Reach et al.|/2005a 2007b). Colored error bars represent photometric data: GALEX far- and near- 
ultraviolet (purple); optical UBV RI or ugriz (blue); near-infrared JHK (green); IRAC 3 — 8um (orange); N-band, IRS 


16 um, and MIPS 24 um (red). IRS 5 — 15m spectra are also shown (solid grey), together with stellar atmospheric 
models (dotted grey), and blackbody approximations to the disk emission (dashed grey). In contrast to the bulk of 
debris disks at main-sequence stars, there is a notable lack of dust significantly cooler than 1000 K. 


dust orbiting polluted white dwarfs, as consistently indicated either by low (or undetected) 24 um 
MIPS fluxes, relative to those at IRAC wavelengths (Figure |4). 

A third shared property among these first half dozen examples is that the infrared excesses 
all become apparent by 2 — 3 um, a trait that is now known to be common (but not universal), 
and hence the most prominent of these are detectable from the ground using K-band spectroscopy 
(2007). Although these first several discoveries are likely 
biased by favoring the bright end of the disk luminosity function, they are highly instructive. The 
emergence of dust emission over the stellar photosphere at these wavelengths indicates there is 
substantial material with T = (2898 K um~!/2) ~ 1000 — 1500 K, demonstrating the circumstel- 
lar debris orbits sufficiently close to pollute the stellar surface. Thus, these findings encouraged 
continued searches for thermal dust emission around metal-enriched white dwarfs 


007A). 


3.2. Theoretical Considerations 


The first several disk discoveries were accompanied by modest theoretical work, but funda- 
mental in the context of Spitzer data, the asteroid disruption model, and the ubiquitous atmospheric 
metals. Optically thin circumstellar dust cannot long survive such close orbits as those inferred 
for blackbody grains such as modeled in Figure [3] The luminosities of typical dusty white dwarfs 
(10-35 SL <107!°L,) are insufficient to remove grains of any size by radiation pressure 
(2008a), and thus in stark contrast to debris orbiting main-sequence stars, there is no 
lower limit to grain sizes in white dwarf disks. Despite this fact, Poynting-Robertson (PR) drag 
will efficiently remove any closely-orbiting, optically thin, small dust grains around these stars 


2006). The timescale for particles of radius a and density p to spiral inward from 
orbital radius r towards a star of luminosity L is 1979) 
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Figure 4: Spitzer IRAC and MIPS photometry of four polluted white dwarfs with circumstellar dust 
[2007b), observed within the first two General Observer cycles (symbols as in Figure[3). As with the earlier cases, the 
infrared excesses begin at 2 — 3 um and share a striking lack of emission from cool dust, as emphasized by the upper 
limits at 24 um for two targets (downward arrows). 
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Here, Qpr is the coupling coefficient for radiation pressure and is ~ | for the case of geometric 
optics (appropriate for micron-sized and larger particles; |Artymowicz|1988). For the luminosity 
range given above, PR drag will typically remove any unshielded, micron-sized and smaller, T ~ 
1000 K dust in less than 10 yr. 

Hence small particles expected to be the most efficient emitters at Spitzer wavelengths must be 
replenished on short timescales or persist in an optically thick configuration. While it is possible in 
principle to replenish an entire, optically thin disk of small particles on decade timescales, because 
at least 20% of white dwarfs are externally polluted over a wide range of metal sinking timescales 
(e.g. Zuckerman et al.2003), this would imply an influx of star-grazing planetesimals, sustained 
over long periods, and is therefore unrealistic (as for the cometary impact model). The tidally- 
disrupted asteroid model circumvents this problem by arranging the debris in a flat and (vertically) 
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opaque configuration 2003). In this geometry, the disk height to width ratio is negligible, 
and only a tiny fraction of particles are subject to the full stellar radiation and drag force. For 
stellar radius R and effective temperature To, dust temperature at radius r is then determined by 


(Chiang & Goldreich} |1997) 
2 \ 1/4 7p\ 3/4 
Tx (= (*) Tete (2) 


It is worth noting that particles in such a ring are cooler than exposed blackbody grains at the same 
orbital radius. 

The flat ring model is fully characterized by stellar R and Tz, the distance to the system d, and 
three free parameters: inner and outer disk radii rjn and out (fixing Tip and Tout), and inclination i 


to the line of sight . The flux from the ring is 2003) 


2Tcosi 
F= 7p 


i: ™ BY(T) ¢ dr (3) 


in 


where By(T) is the Planck function. Substituting x = hv/kT with dust temperature given as in 
Equation [2] the emission from the flat disk can be calculated using 


F120 


dx (4) 


1/3 R2cosi ( 2kTotg Bie hv? [ 3/3 
d? 3hv c2 Jy = er—1 


in 


Figure [5] displays examples of the flat ring model as applied to Spitzer IRAC and MIPS pho- 
tometry for two polluted white dwarfs (2009). As with Figures [3] and [4] these two 
infrared excesses arise between 2 and 3 um, and exhibit the characteristic decreasing flux towards 
24 um. In both cases the disk models reproduce the measured flux excesses rather well, with inner 
disk temperatures similar to that expected for blackbody grains, but with modestly cooler exten- 
sions that better account for the 8 and 24 um emissions. These two stars provide a good showcase 
for the flat ring model, as their silicate emission features (see §3.3 and Figure [6) do not signifi- 
cantly alter their mid-infrared photometry (as for GD 362 in Figure[3). More sophisticated models 
are needed to simultaneously fit infrared spectral features and broad-band fluxes, but the flat ring 
model is the basis. 

The flat disk models naturally possess a modest degree of uncertainty; Tin, Tout, and i are 
somewhat degenerate in how they combine to determine the viewed disk solid angle. Generally 
speaking, the inner disk edges are fairly well constrained by their 2.2 and 3.6 um emission, with 
typical allowed ranges of 10% in Tj, (and hence inner disk radius), and up to 20% at the extremes of 
what might be reasonably fitted (Jura et al.|{2007a). Longer wavelength disk fluxes, however, can 
often be reproduced by either a relatively wide temperature and radii range seen at high inclination, 
or a relatively narrow temperature and radii range seen at modest inclination 
Bergfors et al.|{2014). Despite this degeneracy in the outer disk properties, a significant population 
of high inclination systems is statistically unlikely and hence modest inclinations and narrow radial 


extents are favored (Rocchetto et al.) |2015). 
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Figure 5: Two examples of flat disk models fitted to Spitzer IRAC and MIPS flux measurements (symbols as in Figure 
B). The left hand panels show the full spectral energy distributions together with stellar atmospheric models 


Hippel et al.} |2007 2009), while the right hand panels show the excess fluxes after subtraction of the 


photospheres, and are overplotted with the disk models. It is noteworthy that one of the stars does not exhibit an 
infrared excess in K-band photometry. Both of these particular disks are well modeled with inclinations near 50°; the 
inferred dust temperatures at the ring inner and outer edges are given in each panel, and implying the disks are entirely 
contained within 1 Ro. 


Although the morphological appearance of the Figure[5]models can be approximated by single 
temperature blackbodies (as in Figure|3) of T © (Tin + Tout) /2, the flat ring has important physical 
context. Excluding the 8 um bandpass that partly overlaps the region of silicate emission, excellent 
fits to the thermal continua of all but one of the stars in Figures 3}5|have been achieved using flat 
ring models (2008a). For all these well-modeled systems, the 
resulting inner and outer disk radii lie in the ranges 0.1 —0.2Re and 0.3 —0.4Ro, respectively. 
These outer disk radii can be compared directly to the stellar Roche limit, inside of which tidal 
disruptions are expected. 

The distance at which a lesser orbiting body of density p is tidally split by a greater body of 
density P can be expressed as 
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1/3 
5~a(7) R (5) 


where is a constant which typically has a value in the range 1 — 2.5, depending on multiple 
factors (1999). The classical value of &% = 2.45 is for the case of a rotating, uniform, 
self-gravitating liquid in a circular orbit; whereas in the rigid sphere approximation appropriate for 
stray asteroids the value becomes & = 1.26. For asteroid densities in the range p = 1.0—3.5 gem~> 
and a typical white dwarf of 0.6M. and 0.013 Ro, this results in 6 = 0.8 — 1.2Ro. Therefore, the 
simple flat disk model predicts the circumstellar material orbiting these polluted white dwarfs is 
entirely contained within their Roche limits, and consistent with debris created via tidally disrupted 
asteroids. In contrast, if the dust were not in a self-shielded configuration, 1000 K blackbody 
emitters would lie beyond 1.5 Reo for numerous polluted stars with Teg > 15000 K, and debris at 
such distances would require a distinct physical production mechanism. 

The flat ring model is also consistent with theoretical expectations for a disk dominated by 
solids, evolving at short orbital periods. Without gas pressure, a mature circumstellar disk orbiting 
a white dwarf should not have significant extent above the mid-plane. Any primordial vertical 
height will be flattened via exchange of angular momentum between particles (i.e. dynamical 
relaxation) over many orbital periods, and these will be less than 3.6 hr for material orbiting within 
1 Re around a 0.6 Mq star. Therefore, in the absence of external forces a circumstellar disk orbiting 
white dwarf is likely to resemble the rings of Saturn. 

These closely-orbiting, particulate disks around white dwarfs are likely similar to the plane- 
tary rings of the Solar System in a few important ways. The first similarity is purely geometric, 
as the disk sizes are comparable; the outermost discrete ring of Saturn, the F-ring orbits at 0.2 Ro 
above the planet (center), while the E rings extends to 0.7 Ro. The second analogy is physical, 
as disrupted, minor planetary bodies are a primary source of ring material orbiting all four giant 
planets (1993), as well as the metal-enriched white dwarfs studied by Spitzer. Third, 
particulate rings evolving in isolation spread slowly due to the inefficient transfer of angular mo- 
mentum via viscous forces among solids, and in the case of planetary rings can be further slowed 
by ring moons. Therefore, similar particulate disks at white dwarfs can be relatively long-lived, 
and their evolution dominated by processes at the inner edge, where both radiation and gas produc- 
tion provide additional means to transfer angular momentum. Fourth and finally, owing to these 
fundamental similarities, the observed diversity of rings systems in the outer Solar System is likely 
to be reflected in debris rings orbiting white dwarfs. 


3.3. Infrared Spectra of Disks 


Infrared spectroscopy with Spitzer IRS was successfully obtained for the bulk of dusty white 
dwarfs found prior to the last year of the cryogenic mission; the data for these eight systems are 
shown in Figure (6| It is worth remarking on the ground-breaking sensitivity of these 5 — 15 um, 
R = 100 spectral observations obtained with a 0.85 m telescope (in space), as the 10 um emission 
for all stars but G29-38 have peak fluxes in the range 0.3 — 1.3 mJy. This can be compared with 
one of the most sensitive mid-infrared observations made from the ground, where GD 362 was 
photometrically detected in the N-band at 1.4+0.3 mJy using an 8m telescope 
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Figure 6: Silicate emission features detected in Spitzer IRS low-resolution observations of eight dusty and metal-rich 
white dwarfs; these represent all successful 5 — 15m observations of dusty white dwarfs with IRS, and all but one 
star known prior to Cycle 4. [PG 1015+161 was also observed but no usable signal was recovered, a likely result of 
the short individual exposure times imposed by a bright nearby source.] Each spectrum has been normalized to the 
average of the 5.7 and 7.9 um IRAC fluxes, and smoothed by 15 pixels (0.8 um); the binned data points are highly 
correlated and structures within the broad silicate feature are probably not real. The detections are modest in most 
cases, and the data below 8 um are not shown as this region is typically dominated by noise. Nonetheless, these data 
clearly show broad features, with red wings extending to 12 um that are typical of minerals associated with planet 
formation (Jura et al.|/2009a). All spectra are shown on the same scale for comparison, with the exception of GD 362 
whose spectacular feature requires a larger range. 


2005). Despite this ground-breaking sensitivity, all spectra but one have low to modest S/N 
(2009a), and most disks were confidently detected only between 8 and 12 um where the peak 
in dust emission is typically 20% to 60% stronger than the 6 — 8m thermal continuum. 

Each dataset clearly exhibits strong 10 um silicate emission with a red wing extending to at 
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least 12 um. The measured features are inconsistent with interstellar silicates (Reach et al.|/2005a), 
but are instead typical of glassy (amorphous) silicate dust grains, specifically olivines, typically 


found in the inner Solar System and in evolved solids associated with planet formation (Lisse 
2008). In contrast, the spectra reveal no evidence for emission from polycyclic aromatic 
hydrocarbons. These (carbon-rich) molecular compounds are found in the infrared spectra of the 


interstellar medium 2003), some circumstellar environments (Sloan et al.| {2005}, and 


in comets (Bockelée-Morvan et al.| |1995), with strong features near 6, 8, and 11 um. Hence, 
the circumstellar dust at white dwarfs appears relatively carbon-deficient and rocky 


2009a), and likely similar to the material of the inner Solar System. These findings from infrared 
spectroscopy of debris in situ are consistent with, and almost certainly mirrored by, the carbon- 


poor compositions established via several metal-enriched white dwarfs 
et al 2012} (Farihi etal] 2013B). 

The only target sufficiently bright to be studied spectroscopically with Spitzer at wavelengths 
longer than 15 um is G29-38. obtained a second, low-resolution IRS spec- 
trum of the prototype dusty white dwarf, this time between 5 and 35 um. Interestingly, the repeat 
observations over 5 — 15um revealed no discernible change in the shape of the continuum or sili- 
cate emission, while the longer wavelength data revealed an additional, weaker silicate feature (or 
combination of features) between 18 and 20 um. 

The opaque disk models discussed in the previous section are insufficient to account for the 
observed silicate emission from all dusty white dwarfs observed with IRS. Specifically, the emis- 
sion cannot be modeled in a disk without a vertical temperature gradient or optically thin region. 
Completely optically thin, dust shell models have been used to reproduce the IRS data for G29-38 
(2005a); such models are admittedly attractive as they permit the co-identification 


of various dust species (e.g. minerals, water ice) with the overall shape of the infrared spectra, 
including both continuum and emission components (2009). However, as discussed 
in detail above, optically thin disks are problematic for at least a few theoretical reasons, and 
the body of evidence on polluted white dwarfs suggests many stars have accreted at least 1077 g 
(2010), which is several orders of magnitude 
larger than the ~ 10'° g of small dust grains required for a completely optically thin disk. These 
micron-size disk particles that efficiently emit at Spitzer wavelengths can easily be dwarfed by the 
mass contained in large bodies, but interior to the Roche limit the largest objects will be km-size 
(Davidsson|[1999). 

Alternatively, the flat and optically thick disk models can be readily extended to produce an 
emission feature, by including either a hotter component on the (vertically) outermost layer, or an 
optically thin (radially) outer region 2007b). Such a model is shown in Figure {7| for 
GD 362, utilizing only olivine grains for both the emission components and thermal continuum. 
A similar model was in fact used successfully for G29-38, requiring silicates only 
2009). A further advantage of using the optically thick and geometrically thin disks as the baseline 
model for infrared excess emission around white dwarfs is that these are naturally associated with 
a small but significant fraction of optically thin material. In fact, demonstrate 
that for the dust properties inferred from the opaque disk models (e.g. particle size and density, 
opacity, disk radius), there should be 2 x 10!” g of associated optically thin material. This com- 
pares remarkably well with the total mass of silicates invoked for G29-38 in a completely optically 
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Figure 7: A complete disk model for the thermal continuum and dust emission features at GD 362. Data point with 
error bars are photometry from IRAC and MIPS, while the magenta line is the IRS low-resolution spectrum. The 
blue line is the total model dust spectrum, which is the sum of emission from three disk components: 1) cyan is an 
innermost, opaque disk; 2) red is a radially intermediate zone with a vertical temperature gradient; 3) green is an 


outermost, optically thin region 2007b). The only dust species used in the model are olivines. 


thin model (Reach et al.) |2005a). 


Thus, it appears the flat and optically thick models for circumstellar disks at white dwarfs 
provide an excellent basis for infrared emission features beyond the thermal continuum. The 
extension(s) necessary to account for the silicate emission are relatively straightforward and well 
understood, but it is unclear whether these emitting regions (i.e. disk halos) reside in the radially 
outermost orbits, or at the largest vertical areas. Global models of white dwarf disks indicate that 
optically thin material in the outermost orbital regions is rapidly removed on PR drag timescales 
(2011), implying this area must somehow be replenished with small grains 
to be the prime source of silicate emission, and thus potentially favoring an origin from the top 
layers of the disk. The infrared excess at G29-38 has been present for at least 35 years 
1983), and its strong silicate emission was likely detected in several, independent photometric 
measurements 26 years ago (see Figure Telesco et al.|/1990). While 
not a strong constraint, together with the universally observed silicate emission, this may support 
the self-sustaining production of small (micron- and sub-micron) grains that would otherwise be 
removed on short timescales. 


3.4. Gaseous Debris in Disks 
Early in the Spitzer era, and prior to the publication of any satellite observations of white dwarf 
disks other than G29-38,|Gansicke et al.|(2006) announced the discoveryP|of strong gaseous emis- 
>The identification of metallic emission lines was accidental; its spectrum was flagged in a search for weak spec- 
troscopic features due to very low mass (stellar or substellar) companions to apparently single white dwarfs. 
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Figure 8: Initial spectroscopic observations of SDSS 1228. Upper panel: The SDSS discovery spectrum is shown 
in black, and exhibits emission lines of Call and Fell, as well as photospheric absorption due to Mgt. The red 
line shows a pure hydrogen atmosphere model, where the only expected lines are the Stark-broadened Balmer series. 
Lower panel: Ca It red triplet emission measured at the William Herschel Telescope (WHT) with the ISIS spectrograph 


(Gansicke et al.||2006). The data are shown as points with (small) error bars, while a dynamical, disk model for the 


emission is shown as a solid line. 


sion from a metal-dominated disk orbiting the white dwarf SDSS J122859.93+104032.9 (SDSS 1228). 
This remarkable object proved to be the prototype of a class of dusty and polluted white dwarfs 
with detectable gaseous debris, and this additional observational data has proven highly valuable 

in the study of disk origins and dynamics. 

Figure [8] displays a full optical spectrum of SDSS 1228, including a high signal-to-noise ob- 
servations of the Ca II triplet. Most remarkable are the double-peaked emission lines, whose sub- 
stantial widths are the result of large radial velocity-shifts due to Keplerian rotation. Such lines 
originate from orbiting, heated gas and are typical of accretion disks in cataclysmic variables 


and active galactic nuclei (Horne & Marsh}/1986). However, unlike conventional accretion disks, 
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there is a distinct lack of emission lines from hydrogen; together with the calcium emission, these 
data indicate the disk is significantly depleted in volatiles yet rich in refractory metals. It is no- 
table that Mg II absorption is apparent in the modest resolution (R ~ 2000) SDSS spectrum of a 
Tepe & 22000 K DA star, as photospheric metal detection typically requires Keck or VLT high- 
resolution data for white dwarfs of similar temperature and composition 2005). 


SDSS 1228 remains among the most highly polluted white dwarfs known, with a calcium abun- 
dance close the solar value (Gansicke et al.| {2012). 

There are a few critical, observational properties that apply to the entire class of (single) white 
dwarfs with metallic emission lines. The most fundamental is that all white dwarfs with Ca II 
emission have detectable infrared excesses due to dusty disks (Brinkworth et al.|/2009}|2012); these 
stars are a subset of dusty white dwarfs and the gas component is physically and fundamentally 
related? To date, all white dwarfs with detectable gas emission from their disks also harbor 
substantial atmospheric metal pollution (Gansicke et al.|[2008). This may be an observational bias 
as detection favors 1) high abundances of photospheric metals, 2) bright infrared excesses (large 
disk masses), and 3) strong emission from metallic gas. 

The gaseous emission from debris disks at white dwarfs provides excellent, independent sup- 
port for the conclusions of the infrared data analyses and the asteroid disruption model. Interest- 
ingly, the emitting gas and dust are spatially coincident in all three cases where there are sufficient 
data to model both components in detail (2010a). The observed distribution of veloc- 
ity shifts within the emission lines places the detected gas in white dwarf 
disks in essentially the same orbital regions as those inferred for the dust (Brinkworth et al.| 2009). 


This is a profound empirical result relying only on Kepler’s third law, and confirms all the material 
orbits well inside the Roche limit, significantly strengthening the model of tidally destroyed parent 
bodies for the disk origins. 

Additionally, the shape of the emission lines indicates a clear maximum velocity for the emit- 
ting gas (and thus a minimum orbital radius) on the order of 500kms~!, and corresponding 
roughly to orbits beyond 10 — 20 white dwarf radii. These detected orbital speeds are distinct 
from the 1000 — 3000kms~!, (circular) Keplerian velocities expected for material within 10 stel- 
lar radii of the surface (and inside the observed disk edges); emission from such material has 
not been detected. This last fact may seem surprising, since gas is expected to be produced via 
sublimation near the inner disk edges (as inferred from infrared detection of the emitting solids), 
spiral inward and eventually fall onto the stellar surface. Hence, the behavior of the emission 
lines clearly indicates the detected gas is not produced from sublimation. Further strengthening 
this conclusion is the fact that gaseous emission from disks has been found for the full range of 
stellar effective temperatures where both metals and infrared detected disks are found, down to 


Tere © 13000 K (Farihi et al.) /2012a). 
Furthermore, it is noteworthy that the Ca II line profiles are often asymmetric (Gansicke et al.| 


2008} 2012a; |Melis et al.| |2012), indicating substantial orbital eccentricities for the 


gaseous debris, and variable in shape and strength on timescales from months to years. Figure 
(9| shows the third disk to be discovered with metallic emission lines, Ton 345, with two epochs 


©There is an unfortunate tendency to describe stars like SDSS 1228 as having ‘gas disks’; this is misleading at best, 
and fails to recognize the full properties of these disk systems. 
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Figure 9: The SDSS discovery spectrum of Ton 345 taken 2004 Dec is shown on the left, and on the right is follow 


up data taken 2008 Jan with ISIS on the WHT (Gansicke et al.||2008). The later data trace the shape of the emission 
lines in the SDSS spectrum, and exhibit a clear change in both strength and shape within 3 years. The thin solid lines 
at the bottom of the panels indicate the positions of strong night sky lines. 


of optical spectroscopy separated by 3 years that show clear variation in both the total flux and 
line asymmetries. The possible implications of this variability will be discussed in later sections 
on disk evolution, but it is clear that observations of gas in white dwarf disks will be a highly 
valuable constraint on their nature. At current count there are seven stars with detected emission 
from the Ca II triplet, photospheric metals, and infrared excesses; but only the circumstellar gas 
around SDSS 1228 and Ton 345 has been modeled in detail. 

For these two systems, [Hartmann et al.|(2011| |2014) constructed non-LTE, axially symmetric 
models of gas using a standard & (accretion) disk prescription, where the gas emission is driven 
by viscous heating. The full disk components are built up from concentric rings, where each indi- 
vidual ring is fully solved for radiative transfer and structure, comprising hydrostatic and radiative 
equilibrium, population of atomic levels, metal-line blanketing, and irradiation by the central ob- 
ject. These models are relatively successful at reproducing the broad morphologies of the Ca II 
triplet line profiles, including asymmetries due to eccentric rings. If accurate, these models may 
constrain the chemical composition of the debris; the lack of observed C II emission just redward 


of the Ca IJ triplet is inconsistent with a chondritic mass fraction of carbon (Hartmann et al.,/2011), 
but consistent with the terrestrial value (Hartmann et al.\/2014). The models also yield a strong 


upper limit to the H mass fraction of 1%, consistent with volatile-depleted parent bodies, and 
finds innermost radii for the gas at modestly wider orbits than those found by previous studies 
(Melis et al.|/2010a), supporting previous conclusions that the gas and dust 
coexist spatially. 

Caution is warranted in the above interpretations, however, as a major weakness of the accre- 
tion disk model are strongly over-predicted emission lines of Call H & K, and additional metal 
lines (e.g. Mg I1) that are not observed. In contrast, Melis et al.](2010a) modeled the circumstellar 
gas as ultraviolet-photoionized regions of metals, or Z II regions analogous to HII regions. The 
numerical predictions of the Z 11 model yield total emergent fluxes from emitting gas that are in 
broad agreement with those observed from the Ca II triplet alone. These models therefore over- 
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come the main difficulty encountered by the viscous disk model. It is noteworthy that all models 
for gas origins and behavior in white dwarf disks predict similar gas temperatures, and therefore 
scale heights of at least 10’ cm (Metzger et al.||2012). This clearly implies the bulk of emitting gas 
and dust are vertically distinct, making their radial coincidence remarkable if the two components 
are physically decoupled. These issues are discussed further in §5 on disk origins and evolution. 


4. Disk and Pollution Frequencies 


This section reviews the results of various searches for both atmospheric metals and circum- 
stellar disks, and will attempt to place the previous section into the proper historical and scientific 
context. Notably, the §3 focussed on the discoveries and detections, but the full statistical pic- 
ture requires a discussion of the non-detections and observational sensitivities. The connection 
between the metals and the circumstellar disks has been established empirically by observations, 
and in this section are discussed physical mechanisms contributing to this well-established link. 


4.1. Accretion Rates and History 

Before discussing observational strategies for detecting signatures of planetary debris orbiting 
white dwarfs, it is critical to introduce the concept of inferred accretion rates. In a series of 
pioneering papers, made the first efforts to understand the presence 
and long-term evolution of atmospheric metals in cool white dwarfs, including both gravitational 
settling and episodic accretion. They demonstrated that in a steady-state balance between accretion 
onto the white dwarf surface, and the downward diffusion of heavy elements from the base of the 
outermost layers that are chemically homogeneous with the photosphere, the mass fraction X of 
any metal can be expressed as 


<M 
Mevz 


X= 


(6) 


Here 7 is the exponential sinking or diffusion timescale for a given metal species, M is the steady- 
state accretion rate of the heavy element onto the white dwarf, and M.yz is mass of the convection 
zone. For a given metal, tT can vary significantly as a function of atmospheric temperature and 
composition (see below), but for a given star the sinking timescales for distinct heavy elements 
differ typically by no more than a factor of 2, and by no more than a factor of 3 at the extreme 
(Paquette et a.)[1986). 

Not all polluted white dwarfs have significant convection zones, but at the time of the [Dupuis| 
study, all but one known case were helium-rich (DBZ or DZ) stars which possess 
large convective regions where atmospheric metals are thoroughly mixed. The exception was 
the prototype DAZ, G74-7 (1983), which is sufficiently cool to have a sizable 
convection zone similar to warmer DBZ counterparts. This formalism has since been generalized 
to stars where radiative energy transport is important (or dominant), and M.yz is taken as the mass 
of the outer stellar layers where any heavy elements are well-mixed; this is defined at the base of 
the convection zone or at Rosseland optical depth 5, whichever is deeper (2009). 

Using spectroscopy of white dwarf atmospheres and comparing with spectral models including 
various metal species, the abundance of any heavy element can be determined relative to the main 


PAL 


atmospheric constituent (hydrogen or helium). For any spectroscopically detected metal, its pho- 
tospheric abundance by number is directly proportional to X within the outer layers, and Equation 
(6|can be re-written to infer the (steady-state) rate of infall of an individual element, or summed to 
obtain the total mass accretion rate 


M,=Y m= ye (7) 
j L 


Use of the above equations requires the assumption of steady-state accretion, meaning the rate of 
accretion and diffusion are equal, and equivalent to a constant, downward mass flux for each heavy 


element at every optical depth (Ginsicke et al.| |2012) 
pXv = const (8) 


Here p and v are the mass density and diffusion velocity for a given metal at a specified optical 
layer, determined by atmospheric models and diffusion calculations respectively, and X is mea- 
sured as before via modeling of the spectroscopic data. 

These last two equations have profound importance for the empirical study of exoplanetary 
system assembly and chemistry. They provide the means to calculate the mass accretion rate for 
a specific element within the infalling debris itself, and the ratio of mass accretion rates for any 
two heavy elements is an indirect determination of their intrinsic, relative abundance within the 
parent body (or bodies) of the debris (2010). Assessment of extrasolar minor planet 
compositions using measurements of metal abundances in white dwarfs, and the implications for 
terrestrial planet formation throughout the Galaxy, is outside the scope of this review but readers 
are referred to\Jura & Young]|(2014). While the subtopic of element to element ratios is peripheral 
to this review, the absolute abundances (and hence masses) of metals contained in both the disks 
and the polluted star are also critical. 

Total mass accretion rates must be estimated via extrapolation from detected major elements, or 
treated as a lower limit for undetected yet important species. By mass, the most prevalent elements 
found in polluted white dwarfs are, to date, the same as those that dominate the terrestrial bodies of 


the inner Solar System: O, Mg, Al, Si, Ca, and Fe (McDonoughj |2000;/Visscher & Fegley}|2013); 


C and S are important only for some chondrites and comets; (Lodders & Fegley} |1998} |Lodders| 
2003). Of the few to several hundred known white dwarfs with atmospheric metals, accretion 


rates for these key elements are only known for a handful of stars studied in detail using ultraviolet 
spectroscopy with the Hubble Space Telescope and high-resolution, optical spectrographs on 8m 
class telescope like Keck and the VLT. Although far ultraviolet studies have demonstrated that 
Si 11 (1260/1265 A) is likely the most sensitive indicator (via absorption line strength) of external 
pollution in stars with Ter 2 17000 K (Koester et al.|[2014), a typical observational outcome is the 
detection of this single element. Nonetheless, silicon alone is likely to be a good indicator of total 
accretion rate, as it is one of the four major, rock-forming elements. 

For the bulk of white dwarfs in the optical, only Ca II (3934/3968 A) is detected owing to the 
strength of this atomic transition at relatively cool stellar Tze (the same reason the K line is the 
strongest absorption feature in the Sun). For this reason, accretion rate calculations have been 


22 


historically tied to calcium, thus allowing a diagnostic of ongoing or recent pollution across the 
(essentially) entire population of white dwarfs with metals. Using the infall rate for calcium, one 
can infer the total mass accretion rate as follows 


: 1 Xcqa M. 
Mw Ca 44 cvz (9) 

A Ta 
where A represents the mass fraction of calcium within the accreted material. The early and 
fundamental work on accretion-diffusion in white dwarfs assumed calcium was accreted together 


with all the elements in solar abundance (from interstellar matter; 1993b} 
& Wilken||2006). However, after the first several Spitzer disk detections, it became clear that the 


total accreted mass should be closer to 1% solar by mass 2007a), representing only 
the heavy elements. Currently, using the handful of stars for which detailed measurements of 
at least O, Mg, Si, Ca and Fe are available, the best extrapolation from calcium to total heavy 


element accretion rate is the bulk terrestrial mass fraction, giving A ~ 0.016 (Allegre et al.||1995} 
2010; 2012b). 


Given the utility of the accretion rate diagnostic, it is important to examine when the steady- 
state equations might be valid, and the implications otherwise. While the diffusion timescales for 
all heavy elements in all white dwarf] are universally orders or magnitude shorter than their cool- 
ing ages, they vary substantially as a function of effective temperature and between hydrogen- and 
helium-rich atmospheres (1986). For 25000K > Ter > 6000 K, typical sinking 
timescales are on the order of days to 10* yr for pure hydrogen atmospheres, and from 10 to 10° yr 
for pure helium atmospheres (2009). If a white dwarf is observed to have photospheric 
metals, the probability that it is in a transient phase of accretion, and not a steady-state, is pro- 
portional to the diffusion timescale divided by the lifetime of the pollution event; Ptrans < T/taisk3 
the probability of being in a steady state is ps; = 1 — Ptrans. Although disk lifetimes are somewhat 
uncertain (see §5.3), broadly speaking the DAZ stars are relatively likely to be in a steady state, 
whereas the DBZ and DZ stars quickly approach probabilities that allow for transient accretion 
phases. 

For an isolated instance of a tidally disrupted asteroid, followed by the formation of a disk 
and its subsequent accretion onto the white dwarf, there will be three distinct phases of pollution 
in the stellar atmosphere. During the first (early or build-up) phase, metals will accumulate in 
the outer layers of the star on timescales shorter than diffusion. In this case, the heavy element 
ratios in the star and circumstellar disk are identical as any significant sinking has yet to occur. 
After the passage of one to a few diffusion timescales, the system will approach a steady-state; in 
this condition the metal-to-metal ratios in the circumstellar disk and stellar atmosphere are related 
by the ratio of their sinking timescales (2009). Finally, when accretion from the disk is 
complete, the photospheric metals will remain detectable for up to several diffusion timescales, 
depending on the observational sensitivity. Any calculated mass fractions for elements within 
the debris itself are then dependent on the time between the end of the accretion epoch and the 
observations (generally unconstrained). 


7Radiative forces can compete with and counteract diffusion for stars with Tez 2 25000 K 2014), 


but this review is restricted to those white dwarfs where atmospheric pollution can be confidently inferred. 
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There remain some uncertainties in the state-of-the-art atmospheric models that supply the 
physical constants for accretion rate calculations. The behavior of white dwarf atmospheres at the 
lowest temperatures introduces some difficulty in re-producing extraordinarily broad and complex 


metal absorption features 2011), and the highest (approaching liquid) densities 
found in helium atmospheres leads to uncertainties in the equation of state (Dufour et al.) /2007). 


The recent and highly successful introduction of 3-dimensional atmospheric models has yet to 
tackle chemical diffusion, and the eventual impact remains currently unknown (Tremblay et al.| 
2015). The postulated existence of thermohaline instabilities in accreting white dwarf systems, 
which would potentially increase the inferred accretion rates by orders of magnitude 
lal.||2013), has been independently examined and found to be unimportant, especially in the steady 
state (Koester||2015). From a purely empirical point of view, it is noteworthy that derived accretion 
rates for all DA white dwarfs fall into a broad but single range; stars with temperatures from 5000 
to 25 000 K, cooling ages of 20 Myr to 2 Gyr, diffusion timescales from days to 10° yr, from purely 
radiative envelopes to deep convection zones (2014). Such consistency is unlikely 


to occur if any major physics was missing from the models. 


4.2. Results of Spitzer Searches 


The first robust searches for circumstellar dust around white dwarfs were carried out with 
Spitzer, using IRAC. During its cryogenic lifetime, the instrument was able to image using two 
distinct cameras simultaneously; one at either 3.6 or 4.5 um, and the other at 5.7 or 7.9m. The 
ground-based, photometric studies of G29-38 and GD 362 provided strong indications that those 
two disks emitted the bulk of their radiation at [RAC wavelengths, making it a promising tool to 
search for other systems in this class. Additionally, while these first two examples (and later a 
few other prominent systems) were discovered from the ground, only one other dusty white dwarf 
has ever been detected from the ground beyond the K band (GALEX 1931; {Melis et al.|/2011). 
Importantly, cryogenic IRAC flux measurements provided a broad and accurate characterization 
of the dust temperatures and disk sizes that was not possible based on ground-based data alone. 

There were a few distinct strategies in IRAC observational searches of white dwarfs, most of 
which exploited the fact that both dust and substellar companions will produce detectable infrared 
excesses for a wide range of system parameters®| Early and notable programs included searches 
for dust and companions orbiting a large, brightness-limited sample of white dwarfs 


(2007), high-mass white dwarfs searched for disks due to possible stellar mergers 
jal.| 2006; |Farihi et al.|{2008b), and surveys for dust orbiting metal-rich white dwarfs 
3007; ura etal} 2007) Farihi etal. 2008), 

The basic outcome of the first wave of IRAC searches was two-fold: 1) the only stars found 
to have infrared excesses are those with photospheric metals (2007), and their 
fluxes are consistent with dust rather than substellar companions (Farihi et al.|/2008b), and 2) many 
white dwarfs presumed to be undergoing steady-state accretion do not exhibit infrared excesses, 


including the DAZ prototype (Debes et al.| 2007). The results of these studies indicated that, for 


8The sensitivity of IRAC to spatially-unresolved, substellar companions via 4.5 ym excess emission is competitive 
with state-of-the-art, extreme adaptive optics; yet no closely-orbiting brown dwarf or planetary companions have been 
discovered in well over 200 systems. 
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Figure 10: Total heavy element accretion rate as a function of effective temperature for 108 white dwarfs with known 
(or estimated) Ca abundance that have been observed with Spitzer IRAC. Hydrogen-rich white dwarfs (n = 49) with 
infrared excesses are plotted as red stars, and as black circles otherwise. Similarly, helium atmosphere white dwarfs 
(n = 59) with infrared excess are plotted as blue stars, and as grey circles otherwise. 


those systems where a steady state is likely (i.e. the DAZ stars), the inferred accretion rate was a 
good indicator of stars with and without detectable infrared excess. It remains the case that white 
dwarfs with M, > 3 x 10° gs! are likely to have disks detectable with IRAC, whereas below this 


value infrared excesses are rare 2007a). In a basic sense, this finding appears physically 
sound if the mass (or surface density for optically thin cases) of a circumstellar disk is related to 


the calculated accretion rate onto the stellar surface 201 La). 

Helium-rich stars with metals (e.g. types DBZ and DZ) were initially not targeted in searches 
for dust disks, and accretion rate calculations for these systems were avoided due to the uncertain, 
requisite assumption of a steady state. But in 2007, the DBZ white dwarf GD 40 was confirmed 
to have infrared excess by Spitzer based on a modest, photometric K-band excess in 2MASS 


data 2007a), and GD 362 was shown to be helium-rich (Zuckerman et al.||2007) even 


though its optical spectrum has a DAZ morphology at most spectral resolutions 
2004). Subsequent disk searches with Spitzer began to include all known types of polluted white 


dwarfs, and attempted to analyze all the stars on the same basis of inferred accretion rate 
(2009). Within a few years, the number of disks orbiting helium-rich white dwarfs more than 
doubled (2010b), and currently there are at least ten similar systems. 

The right-hand side of Equation [7| has a numerator that is simply the mass of a given metal 
in the outer layers of the white dwarf, as determined by observations of spectral line strength; 
the steady-state accretion rate is this mass divided its rate of sinking. This formulation yields 
an instantaneous, ongoing rate for those stars actually in accretion-diffusion equilibrium, while 
for stars that are unlikely (or uncertain) to be in a steady state, the equation provides a historical, 


average rate over a single sinking timescale (Farihi et al.}!2012b). In the simplistic case of constant 
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accretion over disk lifetime or diffusion timescale (whichever is shorter), both the instantaneous 
and historical rates among the DAZ and DBZ stars, respectively, should be comparable in the early 
or steady-state phases 2014). 

Figure[10|depicts infrared excess detections and non-detections, as a function of instantaneous 
or average accretion rate, for all the polluted white dwarfs observed with Spitzer during its first 
seven cycles (2014). This sample of 108 metal-enriched stars includes the bulk 
of objects known prior to or independent from those found in the Sloan Digital Sky Survey, and 
represents nearly all polluted white dwarfs specifically targeted in disk surveys 
2014). The figure makes evident a few remarkable properties 


of this sample, and of polluted white dwarfs in general, some of which will be discussed further in 
§5. 
First, the majority of infrared excess detections occur for stars around 10 000 K or warmer, and 


cooling ages less than around 0.5 Gyr (Farihi et al.|/2009). Moreover, among the targets there is a 
decreasing fraction of detected disks with increasing cooling age, and only two possible outliers 


clearly below 10000 K — where the majority of objects are located (2014). This 
result is better visualized as a histogram plotted in Figure [11 |for all targets as a function of effective 
temperature, together with the number in each bin that exhibit an infrared excess. At the warmest 
end, where cooling ages are less than around 50 Myr, all known stars have detectable disks. This 
fraction decreases to the range 0.3 — 0.5 for post-main sequence ages of 50 — 500 Myr, and then 
rapidly moves toward zero for cooling ages approaching 1 Gyr and beyond. 

Second, the sample of known polluted white dwarfs contains a strong metal-detection bias as 
a function of effective temperature. That is, there are few Spitzer targets in Figure [10| with low to 
modest accretion rates and relatively warm temperatures, yet numerous examples at the cool end. 
This property is shared by the entire population of stars whose photospheric metals are detected 
from the ground, and is related to the atmospheric opacity and the ionization balance for the Ca II 
H & K doublet, both of which are a strong function of Ter (Koester & Wilken| {2006). In §2 it was 
mentioned that helium-rich atmospheres are relatively more transparent than their hydrogen-rich 
counterparts, and together with their longer metal residence timescales, resulted in the first dozen 
or so polluted white dwarfs being DZ or DBZ (Sion et al.}|1990b). But it is also the case that metals 
are more readily detected at cooler white dwarfs in general, and only the most dramatic cases are 
detectable from the ground above 20 000 K 2005). 

Therefore, this detection bias favoring extreme cases of pollution among the warmest white 
dwarfs (where radiative levitation is unimportant), combined with the empirical correlation be- 
tween accretion rate and frequency of infrared excess, results in a high fraction of warmer stars 
with disk detections. This conclusion also partly accounts for the dearth of infrared excesses 
among polluted stars with Tes < 10000 K, as the ongoing accretion rates inferred for this subset 
are generally low to modest. 

Third, in Figure[10|the sub-samples of helium- and hydrogen-rich white dwarfs overlap signif- 
icantly in both parameters, yet the former clearly extend to higher accretion rates not represented 
by the latter. There are also a sizable number of DBZ stars with high inferred accretion rates that 


do not reveal infrared excesses, of which there are few among the DAZ types (Girven et al.,/2012). 
These trends indicate that the short diffusion timescales predicted for DA stars are corroborated 
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Figure 11: Histogram of the 108 metal-rich, white dwarf Spitzer targets plotted in Figure The total number 
of objects in each effective temperature bin are shown by unfilled rectangles, while those with infrared excess are 
shown as filled, grey rectangles. A few representative cooling ages are marked along the top of the figure. While 
biases in the detection of atmospheric metals likely inflate the apparent 100% fraction of infrared excesses around 
the warmest white dwarfs, there is a clear decrease in the fraction of disk detections for older post-main sequence 
planetary systems. 


by the data, where detectable disks and high, ongoing accretion rates are strongly correlated. And 
based on the weaker correlation for DBZ white dwarfs, their time-averaged accretion rates are 
likely to reflect their actual, ongoing rates in only 1/3 to 1/2 of cases; for the remaining stars the 
current accretion rates are likely to be significantly lower. 

Thus the statistical picture from the above Spitzer sample is somewhat distorted by target 
selection bias, actual accretion rate, and the capability of IRAC photometry to detect circumstellar 
dust. Regarding sensitivity, the relevant quantity for a detection of spatially-unresolved infrared 
excess 1s the absolute calibration accuracy, including instrumental limitations and uncertainties in 
photospheric models (2008). The IRAC instrument team report 16 errors of 2% in each of 
the four photometric channels (Reach et al.}/2005b); however, photometric observations of isolated 
white dwarfs with Ter > 7000 K (and Rayleigh-Jeans emission at IRAC wavelengths), results 
in deviations from model extrapolations (2007), and within flux ratios 


between channels (Farihi et al.\/2008a), both of which vary by at least 5%. Imperfect extrapolation 
of stellar models to longer wavelengths is also a significant source of uncertainty for infrared 


excess searches, as many (pre-SDSS) white dwarfs lack reliable photometry at any wavelength 
(2009). For these reasons most IRAC fluxes for white dwarfs are 
reported with a 5% calibration error, which is then the limiting factor in any disk search with 
sufficiently high signal-to-noise and accurate short-wavelength fluxes. 

The most subtle excesses confidently identified at white dwarfs, using IRAC are in the range 


30%-45% above the predicted photosphere at both 3.6 and 4.5 um (Farihi et al.) 2009} |2010b} 
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Figure 12: Fractional dust luminosities (t = jp /L,.) for all known white dwarfs with infrared-detected disks 
|chetto et al.|{2015). The central panel shows 7 as a function of cooling age for log g = 8; the dashed and dotted lines 
trace the predicted values for face-on, opaque flat disks with radial extents given by assuming dust grains persist up to 
1500 and 1800 K, and extend out to 1.0 and 1.2 Ro, respectively. At the top is the histogram of detections as a function 
of cooling age, while the center-right panel indicates these as a function of t. The bottom panel plots the ratio of T to 
a predicted maximum (Tgat) for Zin = 1500 K, rout = 1.0RQ. The light grey symbols represent tooo) < 200 Myr, while 
the dark grey symbols correspond to 200 Myr < tego, < 700 Myr. At bottom-right are cumulative histograms of T/Tyat 
for the light and dark grey target subsets, and the distributions for random inclinations are shown as dotted and dashed 
lines, respectively. 


Bergfors et al.) |2014} 2015). This can be contrasted with a strong infrared excess 


like that observed at G29-38, where the disk fluxes are a factor of 5 — 10 stronger than the stellar 
component at the same two IRAC channels. Also in some contrast to the norm, of the 35 disks 
detected by Spitzer to date, there are three cases where an excess is present, but only at one or 
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more of the longest available IRAC wavelengths: G166-58 exhibits excess at 5.7 and 7.9 um only 


(2008a), PG 1225-079 only at 7.9 um, and HS 2132+0941 only at 4.5 um (the 5.7 
and 7.9 um IRAC detectors are not functional post-cryogen; |Bergfors et al.[2014). Thus, while at 
most 50% of infrared-detected disks reveal themselves as excess in K-band photometry or spectra 
(Farihi et al.}/2010b), a further 5%-10% are also missed at 3.6 or 4.5 um, implying a small number 
of disks remain undetectable with warm Spitzer in cycles 6 and beyond. 

Interestingly, it appears that disks with the lowest fractional infrared luminosity are the most 
numerous. This should be expected if the disk population is represented by a normal distribution 
in terms of intrinsic brightness, as random inclinations will not drastically alter the observed pop- 
ulation. Figure [12] demonstrates this result and more, summarizing eight cycles of Spitzer disk 
observations (2015). The central quantity plotted is the fractional infrared (or 
disk) luminosity, and it is shown for all 35 known dusty white dwarfs as a function of cooling 
age,with adjacent number histograms of detection in each parameter. The top histogram reiterates 
the result shown in Figure[I 1]and discussed earlier (in terms of cooling age rather than Tes), while 
the central-right panel demonstrates the bulk of detected disks are at the faint end of t, with the 
number increasing strongly towards the limit. The two panels at the bottom of the figure show that 
1) in contrast to observational bias, disks are preferentially fainter around younger, more luminous 
white dwarfs, and 2) the number of disk detections as a function of T cannot be the result of incli- 
nation (at high significance for the younger population). Together, the frequency and brightness of 
circumstellar disks around white dwarfs suggest a larger population that remains as yet undetected 


(Bergfors et al.}/2014}|Rocchetto et al.,/2015). 


4.3. Independent Studies and Infrared Statistics 


While Spitzer studies have transformed the picture of metals, accretion, and disks around white 
dwarfs, the sample biases are severe, with numerous targets already known to have strong metal 
pollution or ground-based infrared excess or both. There was a single exception during the cryo- 
genic lifetime of IRAC, where data at 4.5 and 7.9um only were acquired for a sample of 122 
white dwarfs (and 2 subdwarfs) based on their K,-band brightness in 2MASS 
(2007). While this sample of stars is unbiased in terms of atmospheric composition, it favors white 
dwarfs with higher effective temperature and lower mass (larger radius), and the resulting sample 
is highly heterogeneous. Regardless, this is the only Spitzer study to probe such a large and diverse 
sample, with two disks identified (among 20 targets known to have atmospheric metals), and thus 
yielding a nominal fraction of 1.6% with detectable disks. 

Considering only metal-rich Spitzer targets, and statistically re-inserting these back into the 
larger, parent samples to which they belong, also provides an estimate of the detectable disk fre- 
quency for white dwarfs in general, resulting in a fraction of 1%-—3% for cooling ages less than 
around 0.5 Gyr, with an uncertain but clear decrease at older ages (2009). If one 
accounts for the youngest cooling ages represented in the Spitzer and parent samples of stars, then 
it is probably safe to say this result applies to cooling ages of 20 to 500 Myr, as neither metal 


pollution nor disk detections have yet to occur for hotter white dwarfs (2014). 

A numerically more robust approach was taken by (2011), by cross-correlating 
16785 DA white dwarfs with Te¢ > 8000 K from the SDSS with the UKIRT Infrared Deep Sky 
Survey (2007). This resulted in the identification of 1884 stars with both H- 
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and K-band data whose ugrizH photometry could be fitted with spectral models to search for 
excess emission in the K band, representing an order of magnitude increase over all white dwarfs 
targeted by Spitzer. The post-main sequence age range in this sample was based purely on the 
fact that an ultraviolet excess is necessary for white dwarf selection within the SDSS. Twelve disk 
candidates were found in total, but follow up observations with Spitzer IRAC revealed that only 
60% of these had infrared fluxes consistent with disks, and more detailed optical spectroscopy 
demonstrated these same stars also had atmospheric metals (2012a). The remaining 
candidates were photometrically confused with background sources in the infrared, or their data 
indicated low-mass companions rather than disks. By assuming that roughly half of all disks are 
not apparent as K-band excess, these studies together suggests that 0.8% of white dwarfs with 
cooling ages less than 1 Gyr have detectable disks. Also noteworthy is that no disk candidates 


were found among several hundred stars with Teg > 25000 K (Farihi et al. {2012a). 


By far the largest and least biased survey of white dwarfs in the infrared was performed by 
WISE, which imaged the entire sky at 3.4, 4.6, 12, and 22 um (Wright et al.|/2010). Auspiciously, 
the first two bandpasses are relatively close matches to the two shortest wavelength IRAC channels, 
and thus direct comparisons between instruments provides verification and eases identification of 
new disks. Despite the all sky capacity of WISE, there is not yet a commensurate catalog of 
white dwarfs (pending GAIA results). Thus the WIRED team opted to cross-correlate the infrared 
database with each of the two largest available white dwarf catalogs: a preliminary version of 


the SDSS DR7 (spectroscopic) white dwarf catalog (Kleinman et al.}/2013), and the older but 
historically invaluable catalog of (McCook & Sion} (1999). 
In the first case, (201 1a) found a detection in at least one WISE band (nearly always 


3.4 um) for 1527 sources associated with a spectroscopically confirmed SDSS white dwarf. Due 
to the nature of the search, just over 1000 of the detections are known or suspected systems with 
low-mass (M dwarf) companions, and where the cool component dominates the infrared emission. 
Therefore, the study first had to determine the achieved efficiency at detecting white dwarfs, and 
to what flux density. The reported 5o sensitivity at 3.4m is 80uJy but 
find the peak in source counts occurs at 50 wy, and that 395 white dwarfs have photospheric 
emission predicted to be above this level. Among their 52 disk candidates, they find six sources 
(three previously known, and three new) above this flux cutoff, and thus a nominal disk fraction of 
1.5%. However, one of these six is a system with an L-type companion (PHL 5038; 
2009), and another was found to have a neighboring (extragalactic) source contaminating the WISE 
photometry (Barber et al.|/2014). Therefore, this sample yields 1% of white dwarfs with detectable 
disks, over a range of cooling ages likely identical to|Girven et al.|(2011), for Tei > 8000 K. While 
this fraction could be somewhat higher, as also cataloged another dozen 


white dwarfs with indeterminate-type excesses and photospheres predicted to be above the flux 
cutoff, the overall survey underscores the need for careful examination and follow up of WISE 
disk candidates, especially below the 5o sensitivity. The relatively large (6”) WISE PSF at 3.4, um 
implies that background contamination is a serious issue, even for relatively bright sources (e.g. 
GALEX 1931;[Debes et al.|2011b} Melis et al [201 1} Rocchetto et al.2015). 

performed archival and ground-based follow up of disk candidates identi- 
fied by|Debes et al.|(201 1a), Of the 52 total candidates, they point out six are previously confirmed 
disk hosts, 17 are problematic due to source blending in SDSS or UKIDSS imaging, and three are 
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misclassified, leaving only half the possibilities. They obtained independent JH imaging for 16 
of the remaining 26 systems, and retrieved archival UKIDSS images for another two (with eight 
stars unexamined). The independent imaging revealed 12 white dwarfs with neighboring sources 
contaminating their W/SE photometry, and should be relatively unsurprising as most targets are 
fainter than the 50 wy (30) statistical flux cutoff defined by the WIRED team. Four stars reveal 
no additional sources in their JH images, and similarly for two others in UKIDSS, leading |Barber| 
to conclude these are all bonafide disks, adding to the six previously known for a 
total of 12. However, they do not re-calculate the number of target photospheres to which WISE 
was sensitive down to the faintest of their newly claimed disks; in fact, this is not possible without 
extrapolation because the WIRED source counts are clearly incomplete below 50 wy 
fal.|/201 1a). Instead [Barber et al.](2014) simply discard 1062 binary suspects and claim an inflated 
fraction of 12 detectable disks in 465 targets, rather than 12 in 1527; the latter is more realistic 
since it is likely many more white dwarf photospheres could have been detected below 30. There- 
fore an accurate fraction of detectable disks will be closer to the latter ratio, and again near 1%. 
This presumes that all of the six new disk claims are real (not guaranteed), but without K-band 
imaging, background galaxies with the reddest colors can be missed. 


The second WISE study (Hoard et al.\!2013) used the significantly brighter source catalog of 
McCook & Sion} (1999), which contains the bulk of all known white dwarfs prior to the SDSS, 


including the nearest populations. Although the catalog contains 2249 entries, after some years of 
vetting by the community (misclassifications, duplications, etc.), 2202 viable targets were iden- 
tified (2007). As discussed previously, reliable flux measurements at shorter wave- 
lengths are necessary to confidently evaluate the presence or absence of infrared excesses, and 
thus [Hoard et al.|(2013) considered only a subset of 1474 white dwarfs detected in 2MASS. After 
a meticulous inspection of numerous individual cases, a total of 12 disk candidates were found, 
seven of which are new discoveries. Here again is a detectable disk fraction that is close to 1%, al- 
though the authors caution that further confirmation is needed for these candidates. While it is pos- 
sible some of these brighter disk candidates are contaminated by additional, unresolved sources, it 
can also be the case that more subtle, authentic disk excesses (e.g. PG 1457—086, HE 0106—3253; 
Farihi et al.2010b) are often missed. However, the[{McCook & Sion](1999) sample is likely to con- 
tain a significant fraction of white dwarfs with Te¢¢ < 8000 K, and based on the combined Spitzer 
studies, will likely result in a smaller measured fraction of detectable disks. 

Most of the aforementioned studies cover a broad range of white dwarf temperatures and cool- 
ing ages, with the exception of the Spitzer studies focused on metal-polluted stars. But based on 
the totality of data, disk detection in the infrared appears more favorable at younger post-main 
sequence ages (2012), but also below temperatures where dust may be rapidly sub- 
limated. There have been two studies that attempted to capitalize on this possible peak in disk 
manifestation over time. 

In cycle 8, Spitzer was used to study a well-characterized sample of 134 relatively young 


white dwarfs with cooling ages in the range 20 — 200 Myr (Rocchetto et al.} 2015). The targets 
were extracted from the Palomar Green (PG) (Liebert et al.||2005) and SPY 2009} 


DA white dwarf samples, where individual stars were selected for only two qualities: 17000 K 
< Topp < 25000 K and F,(1300A) $3x10-" ergs cm~?s~! A~! (this far-ultraviolet brightness 
constraint ensured good signal-to-noise in a parallel, cycle 18 Hubble COS snapshot). Five stars 


a1 


Table 1: Lower Limits for PG and SPY DA White Dwarfs with Infrared Excess 


Ta (K) [Liebert eval ](2005) [Koester etal 2009} 


< 9500 1:16 (6.3%) 0:60 (< 0.02%) 
9500 — 17000 3:88 (3.4%) 9:246 (3.7%) 
17000 — 22500 3:80 (3.8%) 3:138 (2.2%) 
17000 — 25.000 5:113 (4.4%) 3:179 (1.7%) 

Total 9:347 (2.6%) 12:626 (1.9%) 


Note. Entries are based on Table 1 in (2005), and Tables 1 and 2 in|Koester et al.|(2009). Where 


Tere differs, the PG value was adopted for consistency. The SPY column excludes 46 DA+dM binaries as 
temperature determinations are not consistently available. 


in the sample were found to display clear infrared excesses, including two previously known to 
have detectable dust in the infrared, resulting in a nominal fraction of 3.7% 
2015). However, it is important to note that this fraction does not include known binaries, but 
these account for no more than 10% of the PG and SPY samples and thus do not significantly alter 
the fraction of young white dwarfs with infrared excess. 

propose a similar frequency based on an independent K-band excess 
search, sampling stars from the PG catalog alone. The study targeted apparently single DA stars 
with 9500 K < Teg < 22500 K, noting that all but one (then) known dusty white dwarfs lie in 
this temperature range. However, the survey utilized disparate means of infrared excess detection, 
none of which were applied to the bulk of the sample. Of the 117 stars in total, they obtained in- 
dependent JH K, photometry for 78, while for the remaining 39 they relied on 2MASS data alone 
(with most unreliably detected at K, > 14 mag; |Hoard et al.[2007); a subset of 41 targets were also 
observed with low-resolution, 0.8 — 2.5 um spectroscopy. From these mixed data, 11 targets with 
suspected K-band excess were observed with warm Spitzer IRAC, confirming one dust disk. In 
addition to this newly identified infrared excess, the sample included four previously known cases, 


yielding an apparent fraction of 4.3% with infrared detectable debris disks (Barber et al.||2012). 
There are at least two statistical inconsistencies in this result (Rocchetto et al.|/2015). First, there 


are actually 145 single DA white dwarfs in the chosen Ter range (Liebert et al.||2005), and it is not 
clear how the 117 star subset was compiled (including four previously known disks). Second and 


more fundamentally, it is well established that roughly half of white dwarfs with infrared excess do 


not reveal themselves at K-band 
(2013); therefore the frequency reported by should be 
nearly 9% and clearly discrepant with other results. Underscoring this fact, the infrared excesses 
of PG 1126+026 and PG 2328+108 are apparent in IRAC photometry (Rocchetto et al.|/2015), but 
were both missed in K-band spectroscopy (Kilic et al.|/2006} [Barber et al.|/2012). 


A statistically consistent look at the entire PG catalog yields more accurate, but still nominally 
high, infrared disk frequencies across several temperature cuts. Table[I|lists and compares firm 
lower limits for infrared disk detections among the 347 DA stars in the PG sample (Liebert et al.| 


2005), with 626 DA white dwarfs in the SPY survey 2009), suggesting that the 
former catalog may be a statistical outlier among the overall white dwarf population. However, 
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the differences are relatively small and within statistical uncertainties. It may be the case that the 
PG catalog has been more thoroughly probed than the SPY sample and a few more disks await 
discovery in the SPY sample. 

These frequency constraints, derived either partially or completely based on the PG sample, 
are somewhat enhanced when compared to the various surveys discussed above that appear to 
converge around 1%. Furthermore, it is difficult to reconcile this based on stellar temperature or 
post-main sequence system age. For example, the much larger SDSS samples of 
and |Debes et al.| (201 1a) are likely rather similar to the PG stars, since both select white 
dwarfs based on ultraviolet excess. One possibility is statistical variation, and the PG catalog 
currently is then relatively rich in detections by (small) chance. Another interpretation is roughly 
half of detectable disks have been missed in most surveys, either due to lack of necessary wave- 
length coverage, sensitivity (and calibration accuracy), or observations themselves. A uniform 
and consistent analysis of both the PG and SPY samples using W/SE data might provide further 
insight. 


4.4. Gas, Dust, and Metals in the Planetary Paradigm 


The previous sections have focused primarily on infrared observations, where a sizable sur- 
face area of emitting solids is necessary for the detection of circumstellar, planetary material. By 
convention, stars with excess infrared emission are said to have a disk, while those systems con- 
sistent with bare photosphere are often described as lacking a disk. However, this delineation 
(partly limited by language) is sufficiently inaccurate as to obscure the actual scientific picture. 
The Solar System has at least two extensive belts of debris, yet its analog is far from detectable 
with state-of-the-art facilities on the ground and in space (Wyatt| 2008). 

In the same way, gas is detectable at some fraction of dusty white dwarfs, but is presumably 
present in all cases where there is sufficient circumstellar material to produce an infrared excess. 
Gaseous emission from the Ca II red triplet has not yet been observed in the absence of detectable 
dust (2012), suggesting that gas production is more efficient in systems with 
greater masses of circumstellar material. While there are seven white dwarfs where both dust and 
gas emission are detected, and roughly three dozen with infrared excesses in total, the potentially 
strong emission lines favor their discovery and gas detection in white dwarf disks is significantly 
less common. Based on the SDSS DR7 spectroscopic sample of 4636 white dwarfs, all of which 
were examined carefully by eye (2007), five were found to have Call triplet 
emission. Hence the fraction of white dwarfs with detectable gas and dust is close to 0.1%, and 
likely representing less than one tenth of the detected dust disk population 

Gas in circumstellar disks can also be detected via absorption along the line of sight for favor- 
able geometries. For example, the well known edge-on disk around B Pictoris reveals both tran- 
sient and stable gas components that provide critical 
data on dynamics and composition of the parent bodies whose debris is observed. Such detections 
might be relevant for those polluted white dwarfs without appreciable infrared excess, especially 
in the hypothetical case of disks that are dominated by gas-phase debris (2008). Although 
the disks at white dwarfs are most analogous to planetary rings, major differences will occur for 
any substantial gas components; the particulate matter can remain vertically flat while gas pressure 
will cause any vapor to extend far from the disk plane, to roughly 1/100 of the white dwarf radius 
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(Metzger et al.||2012). Importantly, interior to some radius where sublimation of solids is rapid, 


all circumstellar disks orbiting white dwarfs should become completely gaseous. It is not clear if 
the infalling gas settles onto the stellar surface equatorially, or along lines of weak (undetectable) 
magnetic fields, but the latter case might enhance chances of detection via absorption. 

There are at least two, and likely three cases where absorbing circumstellar gas is implied 
via spectroscopic observations of polluted white dwarfs. Perhaps the most intriguing case is 
EC 11246-2923, where the co-addition of nearly three dozen individual, high-resolution spectra 
taken on 11 separate nights reveals two sets of Call H & K absorption; a pair of strong pho- 
tospheric features, and a set of weaker lines (Figure [13} [Debes et al.|2012a). Remarkably, the 
secondary features are blue-shifted from the photospheric velocity by an amount commensurate 
with the gravitational well depth of the white dwarf, and are therefore consistent with circumstellar 
absorption. This particular star does not exhibit an infrared excess in cryogenic IRAC observations 
out to 7.9 um, and thus if the absorption is confirmed as circumstellar (and not interstellar), then 
there are implications for disk frequency and detection. 

Another two examples are for stars with strong infrared excesses, and hence the association of 
additional absorption components with a circumstellar disk is relatively confident. PG 0843+516 
displays overly strong lines of SiIV relative to a photospheric model that matches dozens of ab- 
sorption components for ten total heavy elements, including multiple lines of both Sill and Sill 
(2012). The fact that extra absorption is only observed for the highest ionization 
state, resonance lines of Si indicates a non-photospheric origin. Moreover, the two lines with ex- 
cessive widths are blue-shifted to a velocity consistent with an orbit only a few stellar radii above 
the white dwarf surface. Two similar but weaker analogs are observed at SDSS 1228, where the 
additional absorptions are shallow with a spread of blue-shifted velocities relative to the photo- 
spheric lines. 

Despite these promising detections of gas in both emission and absorption, it is unlikely to 
be detected as readily as long wavelength emission from dust grains. Improved sensitivity in the 
mid-infrared will be realized with JWST, but for imaging studies the only improvements will be 
enhanced spatial resolution (and thus a substantial reduction in photometric confusion), and the 
ability to examine significantly fainter targets. These two potential gains are unlikely to improve 
the scientific picture of disk properties and frequencies, as both Spitzer and WISE are capable 
of targeting sufficient numbers. However, the large mirror of JWST should enable mid-infrared 
spectra of polluted white dwarfs such that their photospheres are detected at sufficient signal-to- 
noise so that more tenuous and subtle excesses might be confidently identified. 

By far the most sensitive indicator of remnant planetary systems at white dwarfs is atmospheric 
metal pollution. In fact, all observational and theoretical work indicates that the presence of a 
circumstellar disk is often (if not always) implied by atmospheric metals, regardless of infrared 
excess detection. One particularly clear demonstration of this fact comes from the combined 
Spitzer IRAC and Hubble COS study of warm DA white dwarfs with metal diffusion timescales 
on the order of days (2015). The two programs observed 
85 targets in common, and find at least 23 stars where accretion must be ongoing to account for 
the presence of atmospheric Sill, but only two of these systems exhibit infrared excesses. The 
unavoidable conclusion is that, over the range of currently detectable metal abundances, roughly 
90% of disks remain undetected in the infrared. 
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Figure 13: An additional absorption component in the Call lines of EC 11246-2923 2012a). The 
secondary features have a velocity shift consistent with a few to several stellar radii above the white dwarf surface, 
consistent with a circumstellar origin. The relatively weak features were detected from the combination of numerous 
individual spectra taken over many nights, and totaling over 9 hours of exposure. 


Another key result of the COS study is the metal pollution frequencies themselves, as this 
ultraviolet search was arguably the most sensitive to date. Assuming a bulk Earth Si abundance, 
at least 23 of the 85 observed white dwarfs must be currently accreting debris at MW, ~ 10°*! gs 
(Koester et al.|[2014). Notably, an additional 25 targets exhibit atmospheric Si that can be sustained 
by radiative levitation at their current Tere, but not within the past 50 Myr at higher temperatures, 
indicating the observed metal is a recent, external contaminant. Thus, the fraction of these stars 
with circumstellar debris is at least 27% and possibly as high as 56%. These results strongly 
corroborate the pioneering studies based on optical, high-resolution spectroscopy, accounting for 


various differences in target properties and instrumentation (Zuckerman et al.||/2003}|Koester et al. 
2005; |Zuckerman et al.}/2010). Based on the more conservative fraction, the COS survey suggests 


that a typical white dwarf with a cooling age of 100 Myr accretes debris for around 30 Myr, and 
at the observed rates would imply total accreted masses of 107!*! g, in excellent agreement with 
studies of older polluted white dwarfs [2010a). Such masses are equivalent to single 
parent bodies with diameters in the range 40-200 km for typical asteroid densities. 

These findings support three important scientific implications: 


1. The fundamental picture of circumstellar debris and pollution at white dwarfs should rec- 
ognize that, in the absence of sustaining mechanisms, where there are atmospheric metals 
there is circumstellar material. In these cases a disk should be considered present and (indi- 
rectly) detected, albeit in a manner typically unavailable for main-sequence stars and debris 
disk in general. Compared to metal pollution, infrared excess is a valuable, but significantly 
less efficient diagnostic of circumstellar matter orbiting white dwarfs, and identifies only the 
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bright end of the disk population (Rocchetto et al.| 2015). 

2. Metals in white dwarfs are likely the most sensitive indicator of planetary debris within the 
inner, terrestrial zones of stars in general. Extrapolating from typical infrared disk detections 
for white dwarfs with known (and relatively high) ongoing metal accretion rates, and as- 
suming disk mass scales with instantaneous accretion rate, those stars with the most modest 
metal pollution (i.e. lowest accretion rates) should be equivalent to disks with Ljp/L, ~ 10~° 
2013a). If broadly correct, this represents two orders of magnitude greater 
sensitivity to warm dust within several AU when compared to similar detections around 
main-sequence stars. 

3. The minimum disk masses implied by the accretion rates and frequencies of polluted white 
dwarfs are equivalent to parent bodies many thousands of times larger than the meter-size 
barrier. Such large planetesimals are predicted as a by-product of major planet construc- 
tion on the main sequence, and the compositions inferred (to date) support terrestrial-like 
analogs. Metal pollution therefore gives the most sensitive lower limit for warm debris con- 


sistent with rocky planet formation 2014), and is comparable to debris disk 


frequencies around young A-type stars that support planet formation processes beyond the 


snow line 2008). 


5. Disk Origin and Evolution 


This section attempts to bring together both observation and theory to reach a relatively com- 
plete picture of the detected and implied circumstellar disks orbiting white dwarfs. Theoretical 
progress in this field has been substantial in the last several years, but the asteroid disruption model 
remains the standard. Because disk properties have been discussed extensively above, this section 
follows a temporal sequence for the creation of disks within white dwarf planetary systems. 


5.1. Planetary System Possibilities 


While it remains to be seen exactly how often exoplanetary systems resemble the Solar Sys- 
tem, it is an excellent starting point for extrasolar asteroid disruptions. The Sun will eventually 
evolve into a white dwarf, after passing through the first ascent (RGB) and asymptotic (AGB) giant 
phases, and these effects have been studied by a few groups. In brief, Mercury will be engulfed on 
the RGB while the orbits of Venus and Earth expand sufficiently due to mass loss that they escape 


direct engulfment on the AGB (Sackmann et al.|/1993); however, tidal forces are likely to entrap 
and destroy Venus but the Earth may remain intact 1996} Schroder et al.|{2008). The 
remnant solar system should be long-term stable (Duncan & Lissauer}|1998). If the main asteroid 


belt is considered as one component of the inner system, then the white dwarf sun will retain 3/5 
of its terrestrial components, 3/4 of the major planets, together with all moons and planetesimal 
belts. 

Fittingly, the Solar System is a place of asteroid destruction. The Kirkwood gaps of the main 
belt are density depletions due to orbital resonances with Jupiter, and objects in these regions 
experience chaotic eccentricity changes until they cross the orbit of Mars or the Earth 
[1995). The body will then gravitationally encounter a planet relatively soon, and 
be ejected, collide with a planet or planetesimal, or impact the Sun (1995). 
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In the latter case, if Sol instead had R © 0.013 Ro typical of a white dwarf, an eccentric asteroid 
would pass inside the Roche limit and be torn apart by gravitational tides. Therefore, the basic 
process that results in polluted and dusty white dwarfs is likely common and implies at least one 
major exoplanet. 

Based on the basic idea of mean-motion (orbital) resonances, there have now been several 
detailed studies examining the feasibility of delivering sufficient planetesimal mass to the surface 
of white dwarfs. was the first to recognize that although the ratio 
of two planetary orbits remains constant during stellar evolution (both expanding so that a = 
a0Mms/Mwa, typically by a factor of 2-4), their critical Hill separation will increase due to the 
reduced mass of the central star. Hence two orbiting bodies that were marginally stable to close 
encounters on the main sequence, can become unstable during the post-main sequence. 
generalized this result for two planetary bodies, demonstrating that instabilities can 
occur for orbits whose separation exceeds the orbit-crossing (Hill criteria) boundary, and for up to 


several Gyr into the white dwarf phase (Mustill et al.||2014). 


The principle of widening orbital resonances or chaotic zones was first applied to a Neptune- 
Kuiper belt analog carried into the post-main sequence (2011). This model predicts 
that a sufficient amount of planetesimal mass 1) remains available in these orbital regions at the 
white dwarf stage and 2) an adequate fraction is scattered into the inner system to account for 
the observed pollution, assuming a range of efficiencies for mass delivered to the stellar surface. 
It is noteworthy that the latter process requires additional planets in the inner system to further 


enhance planetesimal eccentricity so that periastron lies inside the Roche limit (Bonsor & Wyatt, 
2012), and the inner few AU region may be effectively cleared during the giant phases (Mustill et 


Nordhaus & Spiegel} 2013). 

The strength of the extrasolar Kuiper belt model lies in the plentiful, empirical and theoretical 
work on main-sequence star debris disks, and which support these planetesimal disks retaining 
sufficient mass over the necessary, relatively long timescales for white dwarf systems 
(2010). Furthermore, the collisionally evolved mass distributions expected for debris 
disks descending from A-type stars have been shown to be consistent with observed atmospheric 
pollution across the entire range of sinking timescales and metal abundances 2014). 
A similar argument might be made for Oort cloud analogs if they exist in sufficient number, with 
the additional benefit that their eccentricities can more easily be increased to 1 
(2012), especially if the white dwarf is born with a velocity kick due to anisotropic mass loss 
(2015). However, it has been demonstrated that the Oort-like cloud mass that can reach the 
Roche limit of the white dwarf is insignificant when compared to the observed metal abundances, 
and hence cannot be a source for polluted white dwarfs 2014b). 

Outer planetesimal belts similar to the Oort or Kuiper belts should be rich in volatiles, which 
are found to be strongly depleted in the debris orbiting and polluting white dwarfs 
2009). Were similar objects often disrupted near a white dwarf, relatively light elements such as 
carbon and nitrogen that are common in volatile compounds such as ices cannot escape the pristine 


photosphere (Gansicke et al.) |2012) by stellar luminosity, rotation, or magnetic field; there is no 


process favoring heavier or more refractory materials, yet these elements are nonetheless absent or 


markedly deficient relative to solar values (Jura & Young} |2014). Moreover, for bodies of radius 


50km or greater, the bulk volatile content of Kuiper belt analogs should be retained during the 
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luminous post-main sequence phases (Jura & Xu} |2010). 


These compelling findings have led researchers to favor extrasolar, terrestrial-zone parent 
bodies as the origin of external pollution. Exploring a Jupiter-main belt model extrapolated to 
the white dwarf phase, find the 2:1 mean-motion resonance is sufficiently 
widened to perturb and tidally disrupt a significant number of large asteroids over cooling ages of 
10-300 Myr. Scaling from the Solar System, they find a typical system requires 107 — 10* times 
the current mass of the asteroid belt so that a sufficient material is eventually accreted by the white 
dwarf. Given that the progenitors of white dwarfs are typically A-type stars, 2—3 times more mas- 
sive and a few Gyr younger than the Sun, this is possible but still potentially extreme. Another key 
study simulated the influence of a 4 AU planet on a wide belt of planetesimals spaced uniformly 
over 10 AU, with a particular focus on the effect of planet eccentricity and mass on the perturbation 
(and tidal destruction) efficiency (2014). These simulations demonstrate that 
eccentric planets around half the mass of Neptune are the most effective producers of disrupted 
asteroids, relative to more massive giant planets, and their effectiveness increases monotonically 
with eccentricity. 

Interestingly, also find that the mass of the planetesimal disk must 
be a few thousand times larger than the asteroid belt, corroborating the findings of 
and highlighting a current dilemma. While extrasolar asteroid belts of the necessary size 
and mass are consistent with observational limits, and hence not ruled out empirically 
2008), owing to collisional evolution, massive disks within 10 AU are not expected to survive 
beyond a few Myr, and should be depleted by orders of magnitude by 100 Myr, well before the 
post-main sequence (2010). There is currently no proposed 
solution to this puzzle: chemistry favors terrestrial-like parent bodies formed interior to a snow 
line, and planetesimal belt evolution favors the icy, outer reaches. 

It may be the case that asteroid-like belts are relatively common, as are extrasolar Kuiper 
belts. From a formation point of view, a gas giant formed near a snow line can easily prevent the 
coagulation of a planet just inside their orbit, creating an asteroid belt in the terrestrial-zone 
(& Livio}|2013). The exoplanetary systems orbiting the A-type stars HR 8799, Vega, and Fomalhaut 
all appear to share this basic architecture with asteroid belt analogs located 
near (and interior to) their respective snow lines, and with (confirmed or inferred) giant planets just 
exterior. Early simulations demonstrate that the exoplanet HR 8799e will often perturb its exo- 


asteroids towards the inner system (Contro et al.||2015), suggesting that a progenitor architecture 
of polluted white dwarf systems, one which is similar to the Solar System, is not unusual. 


5.2. Formation and Structure of Debris within the Roche Limit 


This section focuses on the parent bodies and their debris as it evolves under the influence 
of tidal disruption and eventual accretion onto the white dwarf surface. A discussion on white 
dwarf disks may default to the debris detected within the innermost region, but planetesimal belts 
that supply the destroyed parent bodies are also debris disks, albeit not yet detected for polluted 
white dwarfs (2014). Several, more distant 
and colder disk candidates have been identified around hot and luminous white dwarfs by 24 um 
excess emission, the most famous example being the central star of the Helix nebula 
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2007). However, their interpretation as genuine debris disks is complicated by binarity and dusty 
outflows associated with their immediate progenitors (Chu et al.|/2011}/Clayton et al.|/2014). 

Tidal destruction for a solid body orbiting a star has not been modeled as fully as have stel- 
lar disruptions around compact objects, but there should be fundamental similarities. The catas- 
trophic event necessitates passage through the Roche radius, which in turn implies e = 0.998 for 
semimajor axes beyond a few AU, and hence the orbit will be virtually parabolic. From energetic 
considerations about half the mass of the fragments will become unbound 
1988), an important outcome that has been overlooked when inferring parent body masses based 
on the observed photospheric metals. Bound material remains on highly eccentric orbits, at least 
initially (1989); rubble pile simulations with only gravitational forces show 
that bound pieces will continue along their original orbits via increasingly filled rings, and that 
some particle clusters require multiple passes to be disrupted 
2014a). 


Assuming only stellar radiation and gravity, it has been shown that rubble piles disrupted into 
collision-less rings will gradually contract under PR drag, on timescales a few to several orders 
of magnitude shorter than the white dwarf cooling age (2015). In this picture, rings 
composed of identically sized particles contract as separate units, with a spread in contraction 
timescales over as many decades as there are particle sizes. This can be orders of magnitude 
longer than a few hundred orbits of the parent asteroid, and can readily exceed a few Myr for cm- 
size particles. Because the known disks have comparable lifetimes based on both theoretical and 


observational studies (Klein et al.|[2010}{Bochkarev & Rafikov}|2011}|Girven et al.|/2012), infrared 


observations are sensitive to any debris orbiting within a few AU but outside the Roche limit. The 
non-detection of contracting rings of debris demonstrates 
that disks form significantly sooner than from PR drag alone, which specifically cannot pull the 
largest fragments (and hence masses) into close orbit. 

Although speculative, it is likely that two additional processes are the key drivers in the forma- 
tion of closely orbiting disks from a highly eccentric stream of debris: collisions and sublimation. 
Micron-sized dust particles known to be present from infrared observations will have their orbits 
altered significantly and rapidly by radiation drag; e.g. 10 um grains with a = 5 AU will have their 
semimajor axis reduced by up to a few percent within a single orbit after disruption 
1979). Therefore at periastron on subsequent passes, there will be a spread of perturbed orbital 
parameters corresponding to the full range of fragment sizes immune to tidal gravity. Keplerian 
velocities within the Roche radius are a few to several hundred kms~!, ensuring that collisions 
will occur at high speed, and likely generate both gas and dust as smaller fragments are produced. 

Gas facilitates the viscous dissipation of energy, and is likely necessary to both compress 
and circularize the debris into the observed compact configurations. Sublimation is critical for 
the eventual accretion of disk material onto the white dwarf surface, but may also be an important 
means of generating gas and viscosity within the first few orbital passes following tidal destruction. 
Because the periastra of disruption events are most likely to occur near the actual Roche radius, 
a typical fragment will spend nearly one hour within this zone (2015). During these 
early stages of disk formation, and prior to significant shielding, myriad fragments and debris will 
be fully exposed to stellar radiation, where they can attain surface temperatures above 1200 K and 


sublimate small grains such as olivines 2011). 
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An interesting possibility is that periastron becomes the focal point for a collisional cascade 
within the first several passages, and then a dominant and stationary site for dust and gas pro- 


duction over at least a few hundred orbits (Jackson et al.||2014). This is precisely where most 


collisions are expected for eccentric orbits (Wyatt et al.| |2010), and where any additional, star- 
grazing bodies would greatly enhance the collisional cross section. In the Solar System, similar 


collisional cascades following disruptions are supported for the rings of Saturn, Uranus, and Nep- 
tune, as evidenced by the size distribution of the largest fragments 2000), and 
moonlets within the Roche limit that are too large to have formed by accretion (Srem¢éevic et al. 
2007). If this simplified picture is basically correct, disks at white dwarfs could form within sev- 
eral tens of orbits of the disrupted parent body, and thus possibly in as little as several hundred to a 
few thousand years. Assuming a Myr total disk lifetime, this relatively short formation timescale 
is consistent with three dozen disk detections in the infrared, where no dust was detected external 
to the Roche radius. 

Major and detailed theoretical efforts have been made to understand the evolution of disks 


orbiting within the tidal radius of white dwarfs 2011} [Rafikov) |201 lalb} [Bochkarev & 
2011} 2012). To begin, it is worthwhile to summarize the possible and 


likely overall disk configuration and constituents. The outer disk should be mainly comprised of 
solid debris, and while particle sizes have little empirical constraint, it is likely that micron- to 
centimeter-sized grains are the norm. Interior to some radius where sublimation is efficient, the 
disk should become entirely gaseous before its eventual accretion onto the white dwarf. Owing to 
the relatively rapid exchange of angular momentum between debris on orbital periods less than a 
few hours, and in the absence of major gravitating bodies (e.g. nearby planets, ring satellites), all 
solid state matter should lie in a geometrically flat plane, with a vertical height comparable to the 
largest disk particles. In contrast, the inner disk will be vertically extended due to gas pressure, 
and hence also rotate slightly slower than the Keplerian speed. 

The particulate disk will have typical orbital velocities of a few to several hundred kms~! 
between the Roche limit and some radius at which grains are effectively sublimated. Despite these 
high velocities, and the fact that collisional timescales are comparable to the short orbital periods 
in optically thick disks, it has been shown that collisions are unlikely to be important in two ways. 
First, the collisional transport of momentum (i.e. viscosity) resulting from differential rotation 
and random motions is orders of magnitude insufficient to account for the inferred rates of disk 


accretion (Farihi et al.|/2008a), for any reasonable particle sizes and disk masses (Metzger et al.| 
2012). Second, energetic collisions will be effectively damped within particulate disks of high 


surface density, and thus cannot be a significant source of gas (Metzger et al.|/2012). 

Disk evolution and structure are strongly influenced by stellar irradiation. The overall shape of 
the dusty component is determined almost exclusively by the radial extent, from near the periastron 
where the parent body is destroyed inward to some radius (r,) where sublimation dominates. Two 
extreme cases that can result in zero radial width, and likely a short-lived and completely gaseous 


disk, are asteroids disrupted either within 1) r<r; < 1Re or 2) r< 1Ro < rs (von Hippel et al. 
2007). Relatively narrow rings should result when these conditions are minimally avoided (see 


Figure(12} 2011; 2015). Furthermore, all dust disks should have sharp outer 


edges, as any optically thin material is quickly dragged into regions of high surface density (Figure 


Bochkarev & Rafikov||2011). Inward transport of mass and momentum is driven by the PR 
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Figure 14: The evolution of an optically thick disk with an initially Gaussian surface density profile under the influence 


of radiation forces (Bochkarev & Rafikov}|2011). The y-axis plots the optical depth as a function of orbital radius, 


represented on the x-axis as integer multiples of the inner dust disk edge. The plot shows the disk profile at four time 
steps measured by the characteristic PR drag timescale T for cm-size particles, with two remarkable features. First, 
an abrupt outer edge develops immediately, and marches gradually inward, re-shaping the disk from the outside even 
while erosion proceeds from the inside. Second, there is a persistent, optically thin tail at the inner disk edge, where 
the PR effect operates and drives inward mass transport. This behavior is generic for all initial disk density profiles. 


effect at the inner edge of the particulate disk (Rafikov) |201 1aj|Bochkarev & Rafikov} |2011), but 
may also be enhanced by the presence of gas (Rafikov,|2011b} |Metzger et al.||2012). Although 


only the innermost solids are exposed to the full starlight, the relatively rapid inward drift and 
sublimation of unshielded particles results in a substantially lower optical depth at this boundary, 
and thus a commensurate increase in the mass of particles drawn inward per unit time 
201 Ta). 

Assuming disks are dominated by dust and rubble, the inward flux of these solids driven by 
radiation pressure is the bottleneck for mass transport. In contrast, the entirely gaseous portion of 
the inner disk should be readily ionized and subject to the magneto-rotational instability (MRD, 
yielding viscous dissipation timescales that are orders of magnitude shorter than those resulting 
from PR drag (Jural[2008} Rafikov|(201 Tb} [Metzger etal. |[2012} Farihi et al 2012). Conservation 
of mass and this dichotomy of inward transport speeds between solid and gaseous debris, implies 
the inner ring of gas will have a surface density orders of magnitude lower than the dust disk. 
This should produce infalling gas that is optically thin, and consistent with observations 


201 1} 2012at |Gansicke et al.) |2012). The final inspiral for disk material may be 
mediated by weak magnetic fields as low as a fraction of a kG (Metzger et al.| /2012), (negligi- 


bly) spinning-up the star and directing the debris stream(s) away from the disk plane onto poles. 
Otherwise the accretion will proceed equatorially, and may induce shocks as the gas orbiting near 
0.01c comes into contact with the white dwarf surface rotating at a few kms~!. Either process 


4] 


may result in short-lived inhomogeneities in the surface distribution of metals (Montgomery et al. 
2008; |Thompson et al.) |2010). 


5.3. Disk Accretion and Lifetime 


The disk evolution models introduced in the previous section provide sets of predictions that 
can be compared to observational data, primarily via metal accretion rates and infrared excesses. 
As a baseline for what follows, consider the models of PR drag on a flat and optically thick disk of 
solids as the baseline evolutionary process. Figure[14|demonstrates that this mechanism erodes the 
inner edge of the dust disk, but actually drives the entire disk structural evolution from the outside 
inward, while simultaneously retaining a clean outer boundary. Therefore, rings of dust should 
narrow as they accrete onto the white dwarf, and the rate of tapering should increase with stellar 
effective temperature and luminosity. Based on the flat disk model, infrared observations appear to 
support these predictions of narrower rings at younger and warmer white dwarfs, although particle 
sizes may also play a role (Figure[12} {Rocchetto et al.[2015). 

Notably, the speed of dissipation onto the central star by PR drag is restricted by geometry 
in the case of opaque rings, and for a given star depends only on the temperature and location 
of the innermost solids; in fact, the accretion rate cannot deviate significantly from a fixed value 
below ~ 10° gs7!; 2011a). In this case, the total disk mass determines its lifetime, 
supplying continuous accretion over at least several Myr for large asteroid masses. For low-mass 
disks that are optically thin, the PR effect operates everywhere simultaneously and for a given 
star the accretion rate depends on the size, density, and orbital radius of particles. Accretion rates 
in such disks are always lower than for their opaque counterparts 2011), 
while disk lifetime is identical to the PR drag timescale, with likely upper bounds approaching 
Myr for cm-size debris at large orbital radii. 

For a range of reasonable assumptions, this baseline mechanism can account for all but the 
highest inferred metal accretion rates. More specifically and importantly, the predictions of PR 
drag on white dwarf disks can account for the full range of instantaneous accretion rates 
(2012), including all the DAZ stars with diffusion timescales less than 1000 yr. 
That is, the PR drag model successfully describes the subset of white dwarfs where ongoing, 
steady-state accretion is likely. Moreover, the disk lifetimes predicted by these models are a good 
match to those inferred from observations based on 1) the fraction of polluted white dwarfs over a 


given range of cooling ages (Zuckerman et al.|/2010; 2014), 2) the mass of metals 


in the outer layers of convective stars divided by typical accretion rates (Girven et al.| |2012), or 
3) the steady-state accretion phase necessary for self-consistent oxygen chemistry (Klein et al.| 


2011). These combined theoretical and empirical data support a picture where the most prominent 
disks have m >> 107° g and hence parent bodies with d >> 40 km (Wyatt et al.||2014). 
Remarkably, the highest inferred accretion rates for polluted white dwarfs — those that are 
inconsistent with mass transport via PR drag — occur only for those stars with helium-dominated 
atmospheres 2012). This fact is epitomized in Figure [15] where a histogram of 
accretion rates for 78 Spitzer targets, calculated assuming a steady state, reveals an upper threshold 
of 2 x 10? gs~! for DAZ stars where ongoing accretion is virtually certain. In contrast the DBZ 
and similar white dwarfs, where metals can persist for 107 = 10° yr, exhibit a clear subset of rates 
beyond the apparent DAZ limit by one to two orders of magnitude. The inescapable conclusion is 
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Figure 15: Histogram of inferred accretion rates for 78 polluted white dwarfs observed with Spitzer (Farihi et al. 
[2012b}, separating the hydrogen (DAZ) and helium (DBZ-type) atmosphere stars. The population of rates exceeding 
M, > 10°? gs! occurs only among the DBZ-type stars where only time-averaged rates can be calculated, and thus 
must be related to their significantly longer metal sinking timescales. 


that some fraction (or possibly all) polluted white dwarfs experience short periods of high intensity 
accretion. Any high-rate accretion episodes must be sufficiently brief that they have eschewed 
detection among hundreds of surveyed DA stars, and also sufficiently extreme that their present- 
day, time-averaged rates still bear the scars (Farihi et al.||2012b). 

There have been a few ideas suggested to produce accretion rates exceeding those resulting 
from PR drag, and all involve gaseous debris. By far the most detailed and developed is the model 
of runaway accretion (Rafikov]|2011b). In this scenario, gas produced by sublimation at the inner 
edge of the particulate disk spreads outward and results in substantial spatial overlap between both 
gas and solids. The sub-Keplerian speed of the gas can impart an additional (aerodynamic) drag 
force on the disk solids, increasing their inward drift speed and leading to an enhanced rate of 
inner gas production via sublimation. Positive feedback results from a rise in the density of gas 
imparting drag on disk solids, which in turn further amplifies the sublimation rate. Successful 
operation of this runaway mechanism is highly dependent on the coupling strength between gas 
and dust, but can produce accretion rates as high as 10!! gs~!, or a few hundred times those driven 
by PR drag alone Metzger et al.) 2012). 

Other methods to produce high-rate accretion require only that a substantial mass of gas is 
produced in a disk system, comparable to the size of a large asteroid. One possibility is an auxiliary 
asteroid and its debris impacting the disk from a previous disruption, where roughly the entire mass 
of the smaller body should become gaseous due to high-velocity collisions between debris streams 
with distinct orbital inclinations (Jura\|2008). A second scenario already discussed is gas produced 
by collisions during the earliest stages of disk formation, on the orbits immediately following the 


destruction of a single asteroid (Farihi et al.| |2012b; |Bear & Soker, |2013). In either case, if 
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sufficient mass in a fully ionized disk of gas results (so that the MRI operates), then standard a 
disk accretion will result, where rates can easily exceed 10!! gs~! by a few 
orders of magnitude. 

In addition to the difference in predicted accretion rates between runaway accretion, and gas- 
dominated accretion model, there is another key distinction. Runaway accretion proceeds initially 
via PR, gradually increasing toward an exponential rise, and ending in a short-lived spike that 
rapidly depletes all the remaining disk, including the bulk of its initial mass (Metzger et al.|/2012). 
Therefore the accretion burst comes at the end of the disk lifetime, and consumes it utterly. In 
contrast, a gas disk produced immediately following a tidal disruption event comes at the beginning 
of the disk lifetime (or middle if a secondary impact), and for at least some fraction of cases, leaves 
behind a dust disk to evolve via PR drag. Figure{I0]shows there are several helium-rich stars with 
super-PR accretion rates (M, > 10° g s_!) that still retain visible dust disks, and thus the data favor 
recent but prior events that allow for persistent debris. It may also be noteworthy that all four 
DAZ stars with detected (overlapping) dust and gas have ongoing accretion rate consistent with 
PR drag. However, if multiple events occur within a single diffusion timescale for the DBZ-type 


stars (Wyatt et al.|/2014), then runaway accretion may operate between planetesimal disruptions. 


5.4. Disk Diversity and Transients 


The above sections present a simplified picture of a single disrupted body of substantial mass 
to account for the majority of polluted white dwarf observations. But in actuality, it is likely that 
white dwarf disks are rather varied, and similar in some respects to the planetary rings of the 
Solar System; such diversity and dynamical activity are in fact empirically supported. Dusty rings 
with large inner holes, gas emission and variability, infrared non-detections and variability, etc. 
all argue for transient phases and active shaping processes. It is often the case that anomalies and 
outliers provide major leaps of scientific understanding. 

Consider the planetary rings of the Solar System, which lie predominantly within the Roche 
limit of their host. They display a stunning variety in size, optical depth, and textures, with moons 
playing fundamental roles as parent bodies, partitions, sculptors, collectors and dispersers of ring 
material (2002). Rings are active environments, where objects are ground down by 
(micro-)meteorite impacts, producing fine dust and plasma (Mitchell et al.||2006), and can harbor 


transient moonlets — caught in a balance between fragmentation and accretion — up to about | km 
in size (2008). Although such detailed phenomena are difficult or impossible to 
witness in the analogous rings orbiting white dwarfs, the processes and structures of planetary 
rings can provide a basic theoretical framework. Beyond any compositional differences between 
their ring constituents, planetary and white dwarf rings are mainly dissimilar in the depth of their 
gravitational wells (i.e. higher velocities and shorter timescales), and stellar irradiation. 

Based on infrared data, the radial extent of the dust rings orbiting white dwarfs is highly 
varied, and not an observational effect induced by low inclination 2015). As 
discussed above, radiative forces sculpt the outer edges of circumstellar dust and drive those radii 
inward, but there are at least three infrared excesses that manifest only beyond 3.6 um (G166-58, 
PG 1225-079, HS 2132+0941 |Farihi et al.[2008a] /2010b} [Bergfors et al.[2014), not to mention that 
only about half of all infrared disk detections become apparent at K band. Because the shortest 
wavelength at which thermal emission arises is strongly determined by dust temperature, these 
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cases correspond to different physical disk properties, and not their radiative environments (which 
scale with Tere). Thus, at least three disks are relatively depleted in dust at inner radii where 
sublimation is unimportant, yet retain higher (presumably optically thick) dust densities further 
out. 

Perhaps stronger evidence of disk diversity is the utter lack of infrared detection for many stars 
where accretion is ongoing 2015). Possible configurations 
that would escape detection include the following, some of which can be applied to detected disks 
with dust depleted inner regions: 1) completely gaseous disks 2008); 2) extremely narrow 
but opaque rings 2010b); 3) optically thin dust (Bochkarev & Rafikov}|2011). In the 


latter two scenarios, the disk (segment) escapes infrared detection due to insufficient surface area 
of the emitting grains, especially in the case of optically thin material (Debes et al.|/201 1b). 
While there is little doubt that disk solids can easily be vaporized by auxiliary impacts due to 
different orbital parameters, it is unclear if multiple events are compatible in both mass and fre- 
quency to remove the bulk of dusty regions that are otherwise detectable. But more concretely the 
ongoing accretion rates determined for warm DAZ stars without infrared excess can be compared 
to predictions for the three above scenarios. Gas disks should dissipate rapidly and result in high 


accretion rates (Metzger et al.|/2012), possibly orders of magnitude larger than observed to date 
(Farihi et al.\/2012b), while optically thin or dense but narrow rings will be accreted at modest rates 
(Bochkarev & Rafikov| 2011) compatible with the non-detections in Figure|10| Furthermore, the 


warm DA stars studied by COS are consistent with a picture where dust often accretes at modest 


rates, and completely gaseous disks dissipate rapidly 2014). Cumulative infrared 


data suggest disks detections become more frequent as their fractional luminosity decreases 
(chetto et al.|/2015), and thus the body of evidence suggests circumstellar dust is likely present but 
often undetectable. 

In contrast, there are at least two cases where spectacularly strong infrared emission cannot 


be explained by flat dust disks alone (GD 56, GD 362; 2007a\b). In these instances, 


warped or flared disk segments have been invoked to increase the emitting surface area, under 
the assumption that they are radiatively induced [2009a). These models account for 
the observations well, but raises the question of why only these two disks suffer from radiative 
warping, but no others that orbit stars of similar luminosity. In analogy with planetary rings, 
it seems more straightforward to cause vertical disk deviations using gravity, and thus requiring 
large fragments equivalent to ring moons. As an example, Daphnis has a diameter of 8 km and 
causes | km vertical disturbances in the A ring of Saturn; this is 100 times larger than the height 
of the unperturbed rings (2009). Somewhat speculatively, these two white dwarfs 
may represent particularly active or young ring systems with persistent, large fragments orbiting 
relatively near to their dusty segments. 

The role of gas is clearly important, yet there is significant uncertainty in its production and 
removal from polluted white dwarf systems. Since all disks that exhibit gaseous emission com- 
ponents are relatively bright in the infrared, it seems evident that gas production is linked to the 
presence of significant solid matter. On the other hand, although sublimation should be occurring 
in all dusty systems and is predicted to spread to the dusty regions (Rafikov||201 1b), the rare nature 
of systems with emission features suggests a particular quality such as total disk mass or youth. 
An impact on a pre-existing disk could enhance the mass of vaporized debris, but calculations 
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suggest that any gas co-orbital with an opaque disk should rapidly condense, becoming essentially 
removed within several hundred orbits (Metzger et al.|/2012). If correct, this implies that sufficient 
gas must be produced on similar timescales in order to be observed. 

All gas emitting disks may display variability in either the strength or shape of their lines, 
but detailed and long-term monitoring data are only just being published 2015). 
SDSS 1617 had emission features that gradually faded over several years after their onset 
(2014), and is consistent with either condensation or o& disk accretion. However, the latter 
should be accompanied by an increased accretion rate (via atmospheric metal abundance) which 
was not observed. Given the significant dusty components orbiting all stars with Ca II emission, 
gas accretion by the MRI may be suppressed due to small dust grains lowering the ionization 
fraction (Metzger et al.| 2012). 

SDSS 1228 has now been monitored for a dozen years and shows intriguing variations in 
the (generally asymmetric) shape of its emission lines, but no significant change in their over- 
all strength. This prototype and brightest example of a debris disk with gaseous emission appears 
to have a quasi-stable, eccentric gas component that is precessing with a period of a few decades 
(2015). Due to the eccentricity of the gas disk, it has been hypothesized the entire 
disk is in the process of circularization and thus continually producing gas via collisions. For 
orbital periods of a few hours at most, this would mean circularization has not occurred within 
20 000-30 000 orbits where there is substantial gas mass and thus viscosity. While speculative, a 
physical mechanism may be required to sustain the eccentric structure, and a large fragment analo- 
gous to a planetary ring moon could have the necessary gravitational influence. Such a large body 
would also be fully exposed above the disk plane, becoming significantly hotter than surrounding 
(shadowed) debris, more prone to sublimation, and thus a possible source for the gas itself. 

Lastly and somewhat surprisingly, at least one dusty white dwarf has been shown to vary 
in the infrared, and interestingly it is also a disk where Call emission is detected. SDSS 0959 
exhibited a drop in its infrared flux by approximately 35% over a period within 300 dy, where the 
flat disk model projects an increase in the size of the dust depleted, inner regions 
2014). Hypotheses include an impact on this inner region and a disk instability, and while both 
are predicted to result in rapid accretion, there has been no change observed yet in either the 
metal absorption or emission features (B. T. Gansicke 2015, private communication); however 
the precise epoch of the infrared decrease may have been prior to the first spectroscopic dataset 
(2012a). Accretion of any vaporized disk mass may be avoided if re-condensation 
occurs after spreading, or if the gas remains stable due to low ionization. While the wide variety 
of disk features and phenomena are not yet fully understood, it is increasingly likely that the disks 
with gas emission are particularly active, and processing recent events that induce variability. 


6. Summary and Outlook 


Empirical and theoretical work support a picture where dynamically active planetary systems 
persist at a significant fraction of all white dwarfs, clearly manifested by atmospheric metals ac- 
creted from circumstellar debris disks. Key points in the development of this field include: 


e Historical firsts were achieved via observations of polluted and dusty white dwarfs. The 
1917 spectrum of vMa 2 is the first evidence collected by astronomers that (with hindsight) 
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indicates the existence of exoplanetary systems. It is also likely that K-band photometry of 
G29-38 represents the first detection of exoplanetary material, circa 1980. 


e Metal pollution describes the presence of atmospheric metals in Ter << 25000K white 
dwarfs. The outer layers of these dense stars should be free of heavy elements, and thus 
require an external — and often continuous — source of infalling metals. Evidence is now 
compelling that circumstellar disks of planetesimal debris are responsible. 


e Infrared observations have unambiguously connected white dwarf metal pollution and de- 
bris disks. Observations and models indicate the disks are flat and opaque, extending from a 
radius where sublimation is efficient to near the Roche limit. Silicate emission features are 
detected in all observed systems, indicating the presence of micron-size grains. 


e Tidal disruption of large asteroids is the standard model for these systems. Debris is not 
detected exterior to the Roche limit, suggesting that disks form rapidly from eccentric parent 
bodies prone to catastrophic fragmentation by gravitational tides. Theoretical modeling 
of the initial disk formation can constrain important parameters such as total masses and 
primordial orbital radii. 


e Planetary systems likely persist at these white dwarfs. The perturbation of large asteroids 
is consistent with at least one major planet, but retaining sufficient planetesimal belt masses 
in the (former) terrestrial zone is an outstanding problem. Kuiper belt analogs can preserve 
ample mass distributions post-main sequence, but to date their volatile-rich chemistry is not 
observed. 


e Disk accretion proceeds primarily through Poynting-Robertson drag. For opaque disks the 
rate is relatively high and fixed by geometry, while for optically thin material, the rates 
are modest; together these account for all systems where a steady-state accretion is likely 
(i.e. warm DAZ stars). The largest inferred accretion rates occur only among time-averaged 
values for DBZ-type stars and indicate high-rate bursts not yet witnessed. 


e Gaseous emission features are rare, but an important probe of dynamical processes and a 
critical real-time view of active planetary systems. The production and behavior of detected 
gas is still somewhat uncertain, but monitoring this material may reveal clues to disk and 
parent body origins, viscous evolution, auxiliary impacts, etc. 


e Disk diversity has been increasing in recent years, and it is now clear that the bulk of 
disks lie below infrared detection thresholds. The sizes and shapes inferred from thermal 
and gaseous emission show compelling evidence for evolution under the influence of white 
dwarf radiation, and possibly gravitational forces from additional bodies. 


White dwarfs are compelling targets for exoplanetary system research, as their polluted sur- 
faces reveal unique and powerful insight into the chemical and physical assembly of exoplanets. 
This is especially true in the case of asteroid analogs, where metal abundances then provide invalu- 
able empirical data for the building blocks of terrestrial exoplanets. Furthermore, white dwarfs are 
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the descendants of intermediate-mass stars that are typically immune to conventional exoplanet 
detection, especially within their inner regions. The review presented here, and specifically the 
research within, provides a broad context to connect origins, formation, and evolution with the 
remnant yet active systems observed today. 

After submission of this review, the first transit events around a white dwarf were reported, 
based on K2 observations of WD 1145+107 (2015). The photometric data are 
fairly compelling, with myriad ground-based follow up observations confirming multiple, broad, 
complex, and evolving transits with periods of approximately five hours and thus near to the Roche 
limit (2016). Furthermore, WD 1145+017 appears highly pol- 
luted with metals and has an infrared excess in WISE, and thus at face value it appears to provide 
striking confirmation of the asteroid disruption model. While it is too early to know what might be 
learned from the continued monitoring of this system, it is likely to provide unprecedented insight 
into disk formation and evolution. Perhaps most interestingly for the bigger picture, it remains to 
be seen if the detection of these transits is most dependent on geometry or timing (or both). 

The near future holds some outstanding prospects for progressing the understanding of disks 
around white dwarfs. Currently there are basically no empirical constraints on the parent bodies 
themselves, apart from their bulk chemical composition. The best understanding of likely source 
regions and sizes comes primarily from theoretical considerations 
2014), and if additional transiting systems are found they may help to independently constrain the 
sizes and orbits of persistent planetesimals. Direct or indirect detection of the parent planetesimal 
belts is possible in the sub-millimeter with ALMA or with future far-infrared space missions, but 
may require a substantial investment to identify favorable systems (Farihi et al.|/2014). In-situ dust 
mineralogy is possible with JWST and should be highly complementary to the bulk abundances 
revealed from atmospheric pollution. Lastly, the detection of giant planets in these systems may be 
possible with JWST, but if any persist within several AU, they will not escape astrometric detection 
by GAIA. 
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